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Abstract 

Transverse momentum spectra and yields of hadrons are measured by the PHENIX collaboration 
in Au + Au collisions at ^/sJ/N' = 130 GeV at the Relativistic Heavy Ion Collider (RHIC). The time- 
of-flight resolution allows identification of pions to transverse momenta of 2 GeV/c and protons 
and antiprotons to 4 GeV/c. The yield of pions rises approximately linearly with the number 
of nucleons participating in the collision, while the number of kaons, protons, and antiprotons 
increases more rapidly. The shape of the momentum distribution changes between peripheral and 
central collisions. Simultaneous analysis of all the pT spectra indicates radial collective expansion, 
consistent with predictions of hydrodynamic models. Hydrodynamic analysis of the spectra shows 
that the expansion velocity increases with collision centrality and collision energy. This expansion 
boosts the particle momenta, causing the yield from soft processes to exceed that for hard to large 
transverse momentum, perhaps as large as 3 GeV/c. 
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I. INTRODUCTION 



Heavy ion reactions at ultrarelativistic energies provide information on strongly inter- 
acting matter under extreme conditions. Lattice QCD and phenomenological predictions 
indicate that at high enough energy density a deconfined state of quarks and gluons, the 
quark-gluon-plasma, is formed. It is expected that conditions in ultrarelativistic heavy ion 
reactions may produce this new state of matter, the study of which is the major goal of the 
experiments at the Relativistic Heavy Ion Collider (RHIC). 

The high energy density state thus created will cool down and expand, undergoing a phase 
transition to "ordinary" hadronic matter. While the tools of choice to study the earliest 
phase of the reactions, and thereby the new state, are probes that do not interact via the 
strong force, such as photons, electrons, or muons, the global properties and dynamics of 
later stages in the system are best studied via hadronic observables. Hadron momentum 
spectra in proton-proton reactions are often separated into two parts, a soft part at low 
transverse momentum (pr), where the shape is roughly exponential in transverse mass rriT = 



yPT + f^oi and a high pt region where the shape more closely resembles a power law. Soft 
production (low px) is attributed to fragmentation of a string between components of 
the struck nucleons, while hard (high pt) hadrons are expected to originate predominantly 
from fragmentation of hard-scattered partons. The transition between these two regimes is 
not sharply defined, but is commonly believed to be near pt ^ 2 GeV/c 

In proton-nucleus (p+A) scattering, these two regimes depend on the colliding system 
size in different ways. The soft production depends on the number of nucleons struck, or 
participating in the collision (Npart)- The number of hard scatterings should increase pro- 
portionally to the number of binary nucleon-nucleon encounters {Ncoii) since these processes 
have a small elementary cross section and may be considered as incoherent. Hard scattering 
also produces color strings which fragment and produce some low px particles, though these 
are much fewer in number than those from the much more frequent soft scatterings. In p+A 
these Npart and Ncoii are connected by a very simple relation, namely Npart = Ncoii + 1- 

In nucleus-nucleus collisions, the number of participant nucleons does not scale simply 
with A, so it is more useful to study scaling with Ncou or Npart- Collisions are sorted 
according to centrality, allowing control of the geometry and determination of Ncou or Npart- 
In heavy ion collisions, one expects secondary collisions of particles (rescattering) to 
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take place, especially among particles with low and intermediate transverse momentum. 
Rescattering may occur among partons early in the collision, and also among hadrons later in 
the collision. Both kinds of rescattering can lead to collective behavior among the particles, 
and the presence of elliptic flow indicates that partonic rescattering is 

important at RHIC. In the extreme, rescattering can lead to thermalization. Rescattering 
has observable consequences on the final hadron momentum spectra, causing them to be 
broadened as shown in this paper. This relates to some of the key questions regarding the 
evolution of the collision: Are the size and lifetime sufficient to attain local equilibrium? Are 
the momentum distributions thermal, and if so, what are the chemical and kinetic freeze- 
out temperatures? Can expansion be described by hydrodynamic models? Momentum 
distributions of hadrons as a function of centrality provide a means to investigate these 
questions and permit extraction of thermodynamic quantities which govern the predicted 
phase transition. 

This paper reports semi-inclusive momentum spectra and yields of vr, K, and p from Au- 
Au collisions at ^/SNN = 130 GeV. The data are measured and analyzed by the PHENIX 
Collaboration in the first year of the physics program at RHIC (Run-1). 

The paper is organized as follows. In Section|n]the PHENIX detectors used in the analysis 
are described. The data reduction techniques using the Time-of-Flight and Drift Chamber 
detectors, along with the corrections applied to the spectra, are described in Section IIIIl 
Functions that describe the shape of the spectra are used to extrapolate the unmeasured por- 
tion in order to determine the total average momentum and particle yield for each particle. 
The overall systematic uncertainties in the spectra are discussed. The resulting minimum 
bias and centrality-selected particle spectra are presented in Section IIVI In Section |V| a de- 
scription of the particle production within a hydrodynamic picture is investigated. For each 
centrality selection, a hydrodynamic parameterization of the rriT distribution is fit simul- 
taneously to the spectra of different species. The data are compared to full hydrodynamic 
calculations. The transition region in px between hard (perturbative QCD) and soft (hy- 
drodynamic behavior) physics is investigated by comparison of extrapolated soft spectra to 
the data. Finally, we study the dependence of the particle yields on the number of nucleons 
participating in the collision. 
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II. EXPERIMENT 



The PHENIX |10|, |ll| experiment at RHIC identifies hadrons over a large momentum 
range, by the addition of excellent time-of-flight capability to the detector suite optimized 
for photons, electrons, and muons. PHENIX has four spectrometer arms, two that are 
positioned about midrapidity (the central arms) and two at more forward rapidities (the 
Muon Arms). A cross-sectional view of the PHENIX detector, transverse to the beamline 
is shown in Figure H Within the two central arm spectrometers, the detectors that were 
instrumented and operational during the ^/snn = 130 GeV run (Run-1) are shown. The 



12j |. The detector systems 



detector systems in PHENIX are discussed in detail elsewhere 
used for the measurements reported in this paper are described in detail in the following 
sections. 



A. CENTRAL ARM DETECTORS 



The central arm spectrometers use a central magnet that produces an approximately 
axially symmetric field that focuses charged particles into the detector acceptance. The two 
central arms are labeled as East and West Arms. The East Arm contains the following 
subsystems used in this analysis: drift chamber (DC), pad chamber (PC), and a Time-of- 
Flight (TOP) wall. The PHENIX hadron acceptance using the TOP system in the East 
Arm is illustrated in Figure |21 where the transverse momentum is plotted as a function of 
the particle rapidity (the phase space) within the central arm acceptance subtending the 
polar angle 6 from 70 to 110 degrees for pions, kaons, and protons. The vertical lines are 
the equivalent pseudorapidity edges, corresponding to \ri\ < 0.35. More details are discussed 
elsewhere 3| • 



1. TRACKING CHAMBERS 



The charged particle tracking chambers include three layers of pad chambers and two drift 
chambers. The chambers are designed to operate in a high particle multiplicity environment. 

The drift chambers are the first tracking detectors that charged particles encounter as 
they travel from the collision vertex through the central arms. Each is 1.8 m in width in the 
beam direction, subtends 90-degrees in azimuthal angle 0, centered at a radius Rdc = 2.2 m, 
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PHENIX Detector - First Year Physics Run 




West Beam View East 

FIG. 1: A cross-sectional view of the PHENIX detector, transverse to the beamhne. Within the 
two central arm spectrometers the detectors that were instrumented and operational during the 
■s/snn = 130 GeV run are shown. 

and is filled with a 50-50 Argon-Ethane gas mixture. It consists of 40 planes of sense wires 
arranged in 80 drift cells placed cylindrically symmetric about the beamline. The wire planes 
are placed in an X-U-V configuration in the following order (moving outward radially): 12 
X planes (XI), 4 U planes (Ul), 4 V planes (VI), 12 X planes (X2), 4 U planes (U2), and 
4 V planes (V2). The U and V planes are tilted by a small ±5" stereo angle to allow for 
full three-dimensional track reconstruction. The field wire design is such that the electron 
drift to each sense wire is only from one side, thus removing most left-right ambiguities 
everywhere except within 2 mm of the sense wire. The wires are divided electrically in the 
middle at the beamline center. The occupancy for a central RHIC Au-|-Au collision is about 
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FIG. 2: The central arm spectrometer acceptance in rapidity and transverse momentum for pions 
(top), kaons (middle), and protons (bottom). 

two hits per wire. 

At the drift chamber location, the field of the central magnet is nearly zero, so the DC 
determines (nearly) straight-line track segments in the r-0 plane. Each track segment is in- 
tersected with a circle at Rdc, where it is characterized by two angles: the angular deflection 
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in the main bend plane, and the azimuthal position in 0. A combinatorial Hough transform 
technique (CHT) is used to identify track segments by searching for location maxima in this 
angular space/citehough. The DCs are calibrated with respect to the event collision time 
measurement (see Section I11B|) . With this calibration, the single-wire resolution in the r-0 
plane is 160 fim. The single-track wire efficiency is 99% and the two-track resolution is 
better than 1.5 mm. 

The drift chambers are used to measure the momentum of charged particles and the 
direction vector for charged particles traversing the spectrometer. The angular deflection 
is inversely proportional to the component of momentum in the bend plane only. Both the 
bend angle and the measured track points are used in the momentum reconstruction and 
track model, which uses a look-up table of the measured central magnet field grid. For this 
data set, the drift chamber momentum resolution is <7p/p = 0.6% © 3.6%p, where the first 
term is multiple scattering up to the drift chambers and the second is the angular resolution 
of the detector. 

In Run-1, there were three pad chambers in PHENIX. Each pad chamber measures a 
three-dimensional space point of a charged track. The pad chambers are pixel-based detec- 
tors with effective readout sizes of 8.45 mm along the beamline by 8.40 mm in the plane 
transverse to the beamline. The first pad chamber layer (PCI) is fixed to the outer edge 
radially of each drift chamber at a radial distance of 2.49 m, while the third layer (PC3) 
is positioned at 4.98 m from the beamline. Both arms include PCI chambers, while only 
the East Arm is instrumented with PCS. The second layer (PC2) is located at an inner 
inscribed radius of 4.19 m in the West Arm and was not installed for Run-1. 

The position resolution of PCI is 1.6 mm along the beam axis and 2.3 mm in the plane 
transverse to the beam axis. The position resolutions of PCS are S.2 mm and 4.8 mm, 
respectively. The PCS is used to reject background from albedo and non-vertex decay 
particles; however, only the PCI is used for the results presented here. The PCI is used 
in the global track reconstruction with the measured vertex position using the beamline 
detectors (see Section IIIB|) to determine the polar angle of each charged track. Both PCI 
and the beamline detectors provide z-coordinate information with a 1.89 mm resolution. 
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2. TIME OF FLIGHT 



The Time-of-Flight detector (TOF) serves as the primary particle identification device 
for charged hadrons by the measurement of their arrival time at the TOF wall 5.1 m from 
the collision vertex. The TOF wall spans 30° in azimuth in the East Arm. It consists of 10 
panels of 96 scintillator slats each with an intrinsic timing resolution better than 100 ps. Each 
slat is oriented along the r-0 direction and provides timing as well as beam-axis position 
information for each particle hit recorded. The slats are viewed by two photomultiplier 
tubes, attached to either end of the scintillator. A ±2(t vr/K separation at momenta up to 
2.0 GeV/c, and a ±2a (7r+K)/proton separation up to 4.0 GeV/c can be achieved. 

For each particle, the time, energy loss in the scintillator, and geometrical position are 
determined. The total time offset is calibrated slat by slat. A particle hit in the scintillator 
is defined by a measured pulse height which is also used to correct the time recorded at each 
end of the slat (slewing correction). After calibration, the average of the times at either end 
of the slat is the measured time for a particle. The azimuthal position is proportional to the 
time difference across the slat and the known velocity of light propagation in the scintillator 
(for Bicron BC404, this is 14 cm/ns). The slat position along the beamline determines 
the longitudinal coordinate position of the particle. The total time of flight is measured 
relative to the Beam-Beam counter initial time (see Section IIIB|) . the measured time in the 
Time-of-Flight detector, and a global time offset from the RHIC clock. Positive pions in the 
momentum range 1.4 < pt < 1.8 GeV/c are used to determine the TOF resolution. The 
timing calibration in this analysis results in a resolution of a = 115 ps.^ 

Particle identification for charged hadrons is performed by combining the information 
from the tracking system with the timing information from the BBC and the TOF. Tracks 
at 1 GeV/c in momentum point to the TOF with a projected resolution aproj of 5 mrad in 
azimuthal angle and 2 cm along the beam axis. Tracks that point to the TOF with less 
than 2.0 aproj were selected. Figure 0] shows the resulting time-of-fiight as a function of the 
reciprocal momentum in minimum-bias Au+Au collisions. 

^ Ultimately, 96 ps results after further calibration, as reported in [l^ . 
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FIG. 3: Scaled Time-of-Flight versus reciprocal momentum in minimum-bias Au+Au collisions at 
■y/sNN = 130 GeV. The distribution demonstrates the particle identification capability using the 
TOF for the Run-1 data taking period. 

B. BEAMLINE DETECTORS 

The beamline detectors determine the collision vertex position along the beam direction, 
and the trigger and timing information for each event. These detectors include the Zero 
Degree Calorimeters (ZDCs), the Beam-Beam Counters (BBC), and the Multiplicity Vertex 
Detector (MVD) and are positioned in PHENIX as shown in Figure lU 

The Zero Degree Calorimeters are small transverse area hadron calorimeters that are 
installed at each of the four RHIC experiments. They measure the fraction of the energy 
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PHENIX Detector - First Year Physics Run 




□ Installed n Active 



FIG. 4: A side view of the PHENIX detector, parallel to the beamline. The beamline detec- 
tors determine the collision vertex position along the beam direction, and the trigger and timing 
information for each event. 

deposited by spectator neutrons from the collisions and serve as an event trigger for each 
RHIC experiment. The ZDCs measure the unbound neutrons in small forward cones {6 <2 
mrad) around each beam axis. Each ZDC is positioned 18 m up and downstream from the 
interaction point along the beam axis. A single ZDC consists of 3 modules each with a depth 
of 2 hadronic interaction lengths and read out by a single PMT. Both time and amplitude 
are digitized for each of the 3 PMTs as well as an analog sum of the PMTs for each ZDC. 

Q 

There are two Beam-Beam counters each positioned 1.4 m from the interaction point, just 
behind the central magnet poles along the beam axis (see Figure The BBC consists of 
two identical sets of counters installed on both sides of the interaction point along the beam. 
Each counter consists of 64 Cherenkov telescopes, arranged radially about the collision axis 
and situated north and south of the MVD. The BBCs measure the fast secondary particles 
produced in each collision at forward angles, with 3.0 < f] < 3.9, and full azimuthal coverage. 
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For both the ZDC and the BBC, the time and vertex position are determined using the 
measured time difference between the north and the south detectors and the known distance 
between the two detectors. The start time (Tq) and the vertex position along the beam axis 
(Zyertex) are calculated as To = (Ti + T2)/2 and Z^ertex = (^i — ^2)/2c, where Ti and T2 are 
the average timing of particles in each counter and c is the speed of light. With an intrinsic 
timing resolution of 150 ps, the ZDC vertex is measured to within 3 cm. In Run-1, the BBC 
timing resolution of 70 ps results in a vertex position resolution of 1.5 cm. 

Event centrality is determined using a correlation measurement between neutral energy 
deposited in the ZDCs and fast particles recorded in the BBCs as shown in Figure El The 
spectator nucleons are unaffected by the interaction and travel at their initial momentum 
from each respective ion. The number of neutrons measured by the ZDC is proportional to 
the number of spectators, while the BBC signal increases with the number of participants. 



III. DATA REDUCTION AND ANALYSIS 



A. DATA REDUCTION 



The PHENIX Level- 1 trigger selected events with hits coincident in both the ZDC and 
BBC detectors, and in time with the RHIC clock. A total of 5M events were recorded 
at ^/sNN = 130 GeV in the ZDCs The collision position along the beam direction 

was required to be within ± 30 cm of the center of PHENIX, using the collision vertex 
reconstructed by the BBC. 

The trigger on both BBC and both ZDC counters includes 92 ± 4% of the total inelastic 
cross section (6.8 ±0.4 barns). A Monte Carlo Glauber model is used with a simulation 
of the BBC and ZDC responses to determine the number of nucleons participating in the 
collisions for the minimum bias events. The Woods-Saxon parameters determined from 
electron scattering experiments are: radius = 6.38±0.06 fm, diffusivity = 0.54±0.01 fm jlG^ . 
and the nucleon-nucleon inelastic cross-section, cr]^+jv = 40±3 mb. An additional systematic 
uncertainty enters the radius parameter since the radial distribution of neutrons in large 
nuclei should be larger than for protons and is not well determined 

The centrality selections used in this paper are 0-5%, 5-15%, 15-30%, 30-60%, and 60-92% 
of the total geometrical cross section, where 0-5% corresponds to the most central collisions. 
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FIG. 5: The event centrality (upper plot) is determined using a correlation measurement of the 

fraction of neutron energy recorded in the ZDCs (vertical scale) and the fractional charge measured 

1 6 

in the BBCs (horizontal scale). The equivalent track multiplicity in each centrality selection is 



Only tracks which are reconstructed in all three dimensions are included in the spectra. 
These tracks are then matched within 2aproj to the measured positions in the TOF detector. 
For each TOF hit, the time, position, and energy loss are measured in the TOF detector. The 
widths of residual distance distributions between projected tracks and TOF hit positions, 
(J proj 7 increase at lower momentum due to multiple scattering. Therefore, a momentum- 
dependent hit association criterion was defined. 

Finally, a requirement on energy loss in the TOF is applied to each track to exclude 
false hits by requiring the energy deposit of at least minimum ionizing particle energy. A (3- 
dependent energy loss cut whose form is a parameterization of the Bethe-Bloch formulajlSl] 
is used, where 

dE/dx ^ (3-^/^ (1) 

and (3 = L/ct, where L is the pathlength of the particle's trajectory from the BBC vertex 
to the TOF detector, t is the particle's time-of- flight, and c is the speed of light. The 
approximate Bethe-Bloch formula is scaled by a factor to fall below the data and thereby 
serve as a cut. The resulting equation is AE = A(3~^^^ where A is a scaling factor equal 
to 1.6 MeV. The energy loss cut reduces low momentum background under the kaon and 
proton mass peaks. The fraction of tracks excluded after the energy loss cut is less than 
5.5%. 

The measured momentum (p), pathlength (L), and time of flight (t) in the spectrometer 
are used to calculate the particle mass, which is used for particle identification: 



2 P 

m = — 



(2) 



The width of the peaks in the mass-squared distribution depend on both the momentum 
and time-of-fiight resolutions. An analytic form for the width of m? as a function of momen- 
tum resolution Op and time of flight resolution is determined using Equation |21 The error in 
the particle's pathlength L results in an effective time width that is included with the TOF 
resolution, cxr, 

a^, = ^m — +4p j^— . (3) 
The momentum resolution of the drift chambers is expressed in the following form 

al=(c,pl]\(C,p'Y, (4) 
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where Ci and C2 are the muhiple scattering and angular resolution terms, respectively. The 
units of S(j)rns are mrad GeV/c. The constant Ki is the momentum kick on the particle from 
the magnetic field and is equal to 87.3 mrad GeV/c. The constant Ci is the width in due 
to the multiple scattering (ms) of a charged particle with materials of the spectrometer up 
to the drift chambers. The C2 term is the angular resolution of the bend angle (a), which is 
the angular deflection in (p of the track segment relative to the radius to the collision vertex. 



Equation m is used in Equation |21 with P = + where uiq is the mass centroid 

of the particle's mass-squared distribution. The mass centroid is close to the rest mass of 
the particle; however due to residual misalignments and timing calibration, the centroid of 
the distribution is a fit parameter in order to avoid cutting into the distribution. The m? 
width for each particle is written as follows: 

Cl ■ Am\l + ^) + C| ■ 4my + C| ■ (^(mg + p2)) 
where the coefficient C3 is related to the combined TOP, 

C, = ^, (8) 

and pathlength contributions to the time width, ax in Equation |H1 From the measured 
drift chamber momentum resolution, Ci = 0.006 and C2 = 0.036 c/GeV. While the TOF 
resolution is 115±5 ps, the pathlength uncertainty introduces a width of ~ 20-40 ps, so 145 
ps is used for o"r in C3. 

The pions, kaons, and protons are identified using the measured peak centroids of the 
distribution and selecting 2a bands; shown as shaded regions in Figure IHl for two different 
momentum slices. The 2a bands for pions and kaons do not overlap up to pt =2 GeV/c. 
The protons are identified up to = 4 GeV/c. By studying variations in the centroid 
and width before the particle identification cut is applied, the uncertainty in the particle 
identification is estimated to be 5% for all particles. 

Kaons are depleted by decays in fiight and geometrical acceptance. For the low mo- 
mentum protons, energy loss and geometrical acceptance cause a drop in the raw yield for 
Pt < 0.5 GeV/c, as seen in Figure |21 
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FIG. 6: The mass-squared distributions of positive pions, kaons, and protons for two different 
momentum slices. The momentum sHce 0.7 < p < 0.8 GeV/c is the upper panel and 1.3 < p < 1.5 
GeV/c is the lower panel. The shaded regions correspond to the 2a particle identification bands 
based on the calculated mass-squared width, the measured mass-squared centroids, and the known 
detector resolutions. 

The remaining background contribution was determined by reflecting the track about the 
midpoint of PHENIX along the beamline and repeating the association and PID cuts used 
in the TOF detector. This random background was evaluated separately for each particle 
type. The background contribution is ~ 30% for the kaon spectra at 0.2 < px < 0.4 GeV/c 
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and defines the low px limit in the spectra. The background is < 5% in all other cases, 
and negligible above 0.8 GeV/c in the measured momentum range in this analysis. The 
background was not subtracted but is instead treated as a systematic uncertainty. This 
uncertainty is 2, 5, and 3% for pions, kaons, and protons, respectively, at < 0.6 GeV/c 
and is negligible at higher momenta. 

B. ANALYSIS 

The raw spectra include inefficiencies from detector acceptance, resolution, particle decays 
in flight and track reconstruction. The baseline efficiencies are determined by simulating 
and reconstructing single hadrons. Multiplicity dependent effects are then evaluated by 
embedding simulated single hadrons into real events and calculating the degradation of the 
reconstruction efficiency. 

1. CORRECTIONS: ACCEPTANCE, DECAYS IN FLIGHT, AND DETECTOR RE- 
SPONSE 

The corrections for the finite detector aperture, pion and kaon decays in flight, and the 
detector response are determined using single particles in the the GEANT [19] simulation 
of the detector. All details of each detector are modeled, including dead channels in the 
drift chambers, pad chambers, and Time-of-Flight detector. All physics processes are au- 
tomatically taken into account, resulting in corrections for multiple scattering, anti-proton 
annihilation, pion and kaon decays in flight, finite geometrical acceptance of the detector, 
and momentum resolution, which affects the spectral shape above 2.5 GeV/c. 

The drift chamber simulated response is tuned to describe the response of the real drift 
chambers on the single-wire level. This is done using a simple geometrical model of the drift 
chamber and the straight-line trajectories of particles from the zero-field data. This simple 
model of the drift cell in the drift chamber is sufficient to describe the observed drift distance 
distribution, the pulse width, the single wire efficiency, and the detector resolution. The TOF 
response is simulated by smearing the true time of flight using a Gaussian distribution with 
a width as measured in the data. 

Figure [7| shows the momentum dependence of the residual distance between projected 
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tracks and TOF hits for the real (sohd hne) and simulated (dashed) events. These residuals 
are parameterized in the azimuthal angle and the beamline direction separately for data 
and simulation. For each case, tracks that fall outside 2(t of the parameterized width are 
rejected, thus allowing use of the Monte Carlo to evaluate the correction for the 2a match 
requirement for real tracks. 




FIG. 7: Comparison of the momentum-dependent residuals of DC tracks matched to TOF hits in 
azimuthal angle ^ (left) and z (right) between data (solid) and simulation (dashed). 

A fiducial cut is made in both the simulation and in the data to ensure the same fiducial 
volume. The systematic uncertainty in the acceptance correction is approximately 5%. 

The simulated distributions arc generated uniformly in p^, <P: and y. For each hadron, 
sufficient Monte Carlo events are generated to obtain the correction factor for every measured 
Pt bin. The statistical errors from the correction factors were smaller than those in the data 
and both are added in quadrature. 

The distribution of the number of particles generated in each pt slice, dN/dp^, is the 
"ideal" input distribution without detector and reconstruction effects. This distribution is 
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normalized to 27r and 1 unit of rapidity. After detector response and track reconstruction, 
the output distribution is the number of particles found in each px slice. The final correc- 
tions are determined after an iterative weighting procedure. First, the flat input and output 
distributions are weighted by exponential functions for all particles using an inverse slope of 
300 MeV. The ratio of input to output distributions is determined as a function of momen- 
tum. In each px slice, the corresponding ratio is applied to the data. The corrected data are 
next fitted with exponentials for kaons and protons (see Equation [TT|) . and a power-law for 
the pions (see Equation 1^1). The original fiat input and output distributions are weighted 
by these resulting functions. The procedure is repeated until the functions remain constant 
in their parameters. The weighted input and output distributions are divided to produce 
acceptance correction factors. The corrections are larger for kaons due to the decays in 
flight. The statistical error in determination of the correction factor is added in quadrature 
to the statistical error in the data. 



2. HIGH TRACK-DENSITY EFFICIENCY CORRECTION 

A final multiplicity dependent correction is determined using simulated single-particles 
embedded into real events. This correction depends on both the quality of the track recon- 
struction in a high multiplicity environment and the type of particle measured. 

Depending on the centrality of the event, the correction factor is determined for each 
particle in the raw transverse momentum distribution and is applied as a weight. The final 
efficiency corrections are shown in Figure |Hl where the correction for pions is shown as 
solid circles and for (anti)protons as open circles. The horizontal axis ranges from the most 
central to the most peripheral events in increments of 5%. The systematic uncertainty in 
the multiplicity efficiency correction is 9%. 

The difference between pions (solid) and (anti)protons (open) is due to the different TOF 
efficiencies for each particle (protons are slower than pions). In a small fraction of cases two 
particles may hit the same TOF slat at different times, and the slower particle is assigned 
an incorrect time. The particle will then fall outside the particle identification cuts. This 
effect depends on the type of particle. 

For each particle, two curves are shown, representing the DC tracking inefficiency for two 
types of tracks: fully reconstructed and partially reconstructed tracks. Fully reconstructed 
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tracks include XI and X2 sections. In a high track-density environment, tracks may be 
partially reconstructed or hits may be incorrectly associated. There are two cases when 
this incorrect hit association occurs. In the first case, the direction vector in the azimuth 
prevents the track from pointing properly to the PCI detector, and the correct hit cannot 
be associated. In the second case, the track is reconstructed properly, but there are two 
possible PCI points. If no UV hits are found, then the wrong PCI point can be associated 
to the track and the track's beamline coordinate is mis-reconstructed. In both of these cases, 
the track fails the matching criteria in the TOP detector and is lost. 
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FIG. 8: The multiplicity dependent efficiency correction for pions (solid) and (anti)protons (open) 
for two types of tracks. The upper set of points correspond to fully reconstructed tracks in the 
drift chambers; while the lower set of points correspond to partially reconstructed tracks in the 
drift chambers. 



3. DETERMINING THE YIELD AND MEAN pr 

The dN/dy and {px) are determined using the data in the measured region and an 
extrapolation to the unmeasured region after integrating a functional form fit to the data. 
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A function describing the spectral shape is fit to the data, with varying pt ranges to control 
systematic uncertainties in the fit parameters. The fitted shape is extrapolated, integrated 
over the unmeasured range, and then combined with the measured data to get the full yield. 
Two different functions are used to estimate upper and lower bounds for each spectrum. The 
average between the upper and lower bounds is used for dN/dy and {pr)- The statistical error 
is determined from the data, and the systematic uncertainty is taken as 1/2 the difference 
between the upper and lower bounds. 



For pions, a power-law in px (Equation ^ and an exponential in mj'{= + '^o) 
(Equation ITIHl are fit to the data. For kaons and (anti)protons, two exponentials, one in px 
(Equation ITT|) and the other in uit are used. The pt exponential provides an upper limit 
for the extrapolated yield, which is most important for the (anti)protons. The power-law 
function has three parameters labeled A, po, and n in Equation El The exponentials have 
two parameters, A and T. 



2npTdpTdy \po + pr 



2'KmTdmTdy 
d'^N 
2npTdpTdy 



Ae-^^^l^ (10) 
Ae-P^l^ (11) 



C. SYSTEMATIC UNCERTAINTIES 



In Table lU the sources of systematic uncertainties in both (p-r) and dN/dy are tabulated. 
The sources of uncertainty include the extrapolation in p^, the background, and the Monte 
Carlo corrections and cuts. The uncertainty in the Monte Carlo corrections is 11% and 
includes: the multiplicity efficiency correction of 9%, the particle identification cut of 5%, 
and the fiducial cuts of 5%. The uncertainties in the correction functions are added in 
quadrature to the statistical error in the data. Background is only relevant for p^ < 0.6 
GeV/c in the spectra. 

The total systematic uncertainty in the (pt) depends on the extrapolation and back- 
ground uncertainties; the uncertainties are 7%, 10%, and 8% for pions, kaons, and protons, 
respectively. The overall uncertainty on dN/dy includes the uncertainties on (pt) in addi- 
tion to the uncertainties from the corrections and cuts; the uncertainties are 13%, 15%, and 
14% for pions, kaons, and protons, respectively 
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The hadron yields and (pt) values include an additional uncertainty arising from the 
fitting function used for extrapolation to the unmeasured region at low and high pt- The 
magnitude of the extrapolation is 30 ± 6% of the spectrum for pions, 40 ± 8% for kaons, 
and 25 ± 7.5% for protons [2^. The systematic uncertainty quoted here is taken as 1/2 the 
difference between the results from the two different functional forms. 

TABLE I: The sources of systematic uncertainties in {px) and dN/dy. 



TT (%) K (%) (anti)p (%) 

Extrapolation 6 8 7.5 

Background {pT < 0.6 GeV/c) 2 5 3 

(pt) total 7 10 8 

Corrections and cuts 11 11 11 

dN/dy total 13 15 14 



The momentum scale is known to better than 2%, and the momentum resolution affects 
the spectra shape, primarily for protons, above 2.5 GeV/c. The momentum resolution is 
corrected by the Monte Carlo. As other sources of uncertainty on the number of particles 
at any given momentum are much larger, momentum resolution effects are neglected in 
determining the overall systematic uncertainty from the data reduction. 

IV. RESULTS 

A. TRANSVERSE MOMENTUM DISTRIBUTIONS 

The invariant yields as a function of pt for identified hadrons are shown in Figure while 
Figure ITUl provides the centrality dependence of the spectra. The spectra are tabulated in 
Appendix El The vr^, K^, p, and p invariant yields for the most central, mid-central, 
and the most peripheral collisions, were reported previously [2^. Pion and (anti-)proton 
invariant yields are comparable for pt >1 GeV in the most central collisions. 

As can be seen already from Figure^] all the spectra seem to be exponential; however, 
upon closer inspection, small deviations from an exponential form are apparent for the more 
peripheral collisions. The spectrum in the most peripheral collisions is noticeably power- 
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FIG. 9: The spectra of positive particles (left) and negative (right) in minimum-bias collisions 
from Au+Au collisions at ^/sNN =130 GeV. The errors include both statistical and systematic 
errors from the corrections. 

law-like when compared to the more exponential-like spectrum in central collisions. This is 
especially apparent for the pions. The effect can be seen more clearly in the ratio of the 
spectra for a given particle species in two different centrality classes. Such ratios are shown 
in Figure im for the 5% central and the most peripheral positive spectra (60-92% centrality). 
The ratios for protons and antiprotons as well as for tt~^ have a maximum at intermediate 
Pt and are lower both at low and high pt- The kaon shape change is not very significant, 
given the current statistics. 

The change in slope at low-pr in central collisions compared to peripheral is consistent 
with a more substantial hydrodynamic, pressure-driven transverse flow existing in central 
collisions, since the increased boost would tend to deplete particles at the lowest px (see 
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23l | . It is in contrast 



Section ITV CI) . This is observed at lower energies at the CERN SPS 

n 

to results obtained at the ISR [2J| for p+p collisions at ^/s = 63 GeV, where a shallow 
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FIG. 10: The hadron spectra for five centralities from the most central 0-5% to the most peripheral 
60-92% at ^/sNN = 130 GeV. Errors include both the statistical and point-by-point error in the 
corrections added in quadrature. 

maximum or minimum exists at low pt (in the range 0.3 — 0.6 GeV/c). 

1. FEED-DOWN CONTRIBUTION TO p AND p FROM INCLUSIVE A and A 

Inclusive A and A transverse momentum distributions have been measured in the west arm 
of the PHENIX spectrometer using the ^tracking detectors (DC, PCI) and a lead-scintillator 
electromagnetic calorimeter (EMCal) j25|. The invariant mass is reconstructed from the 
weak decays A p -\- and A —>-p + tt"*". 

The tracks from the tracking detectors are required to fall within 3cr of EMCal measured 
space-points. The EMCal timing resolution of the daughter particles is ~ 700 ps. Using the 
DC momentum and the EMCal time-of-fiight, the particle mass is calculated, and protons, 
antiprotons, and pions are identified using 2a momentum-dependent mass-squared cuts. A 
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FIG. 11: The ratio of the most central to the most peripheral yields as a function of pT for pions 
(top), kaons (middle) and (anti)protons (bottom). 

clean particle separation is obtained using an upper momentum cut of 0.6 GeV/c and 1.4 
GeV/c for pions and protons, respectively. The momentum is determined assuming the 
primary decay vertex is positioned at the event vertex and results in a momentum shift of 
1-2% based on a Monte Carlo study. 

Using all combinations of pions and protons, the invariant mass is determined. The 
mass distribution shows a A peak on top of a random combinatorial background, which is 
determined by combining protons and pions from different collisions with the same centrality. 
A signal-to-background ratio of 1/2 is obtained after applying a decay kinematic cut on 
the daughter particles. Fitting a Gaussian function to the mass distribution in the range 
1.05 < mpTT < 1.20 GeV/c^, 12000 A and 9000 A are observed, with mass resolution 5m/ m ^ 
2%. The reconstructed A and A spectra are corrected for the acceptance, pion decay-in- 
fiight, momentum resolution, and reconstruction efficiency The systematic uncertainty 
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on the pt spectra is 13% from the corrections and 3% from the combinatorial background 
subtraction. The feed-down contributions from heavier hyperons and Q are not measured 
but are estimated to be < 5%. 

In Figure^! the transverse momentum spectra of inclusive protons (left) and antiprotons 
(right) are shown with the inclusive A and A transverse momentum distributions. The solid 
points are the (anti)proton spectra after the feed-down correction from A and A weak decays. 
From here forward, the data that are presented and discussed are not corrected for this feed- 
down effect; inc usive p and p yields are given. More details on the A and A measurement 
are included in 2^ . 
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FIG. 12: For minimum-bias collisions, the inclusive A, inclusive p, and direct p transverse momen- 
tum distributions are plotted together in the left panel. The equivalent comparison for inclusive 
A, p, and direct p transverse momentum distributions is in the right panel. 



B. YIELD AND (pr) 

The yield, dN/dy, and the average transverse momentum, (pr), are determined for each 
particle as described in the preceding section and have been previously published in 21 1. 
For each centrality, the rapidity density dN/dy and average transverse momentum (pj-) are 
tabulated in Tables HH and IIIH respectively. 
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The Npart and Ncoii in each centrahty selection are determined using a Glauber-model 
calculation in 2^ . The resulting values of Npart and Ncoii are also tabulated in Table lUl 
(See Appendix ^for more detail). The errors on Npart and Ncou include the uncertainties 
in the model parameters as well as in the fraction of the total geometrical cross section 
(92% ±4%) seen by the interaction trigger. The error due to model uncertainties is 2% |26|. 
An additional 3.5% error results from time dependencies in the centrality selection over the 
large data sample. 

TABLE II: The dN/dy at midrapidity for hadrons produced at midrapidity in each centrality 

class. The errors are statistical only. The systematic errors are 13%, 15% and 14% for pions, 

kaons, and (anti) protons, respectively. The Npart and Ncou in each centrality selection are from a 

I I 

Glauber-model calculation in 26], also shown with systematic errors based on a 92ib4% coverage. 



Npart 
Ncoll 


0-5% 

347.7 ± 10 
1008.8 


5-15% 
271.3 ±8.4 
712.2 


15-30% 
180.2 ±6.6 
405.5 


30-60% 
78.5 ±4.6 
131.5 


60-92% 
14.3 ±3.3 
14.2 




276 ±3 


216 ±2 


141 ± 1.5 


57.0 ± 0.6 


9.6 ±0.2 


7r~ 


270 ± 3.5 


200 ±2.2 


129 ± 1.4 


53.3 ±0.6 


8.6 ±0.2 


K+ 


46.7 ± 1.5 


35 ± 1.3 


22.2 ±0.8 


8.3 ±0.3 


0.97 ±0.11 


K- 


40.5 ±2.3 


30.4 ± 1.4 


15.5 ±0.7 


6.2 ±0.3 


0.98 ±0.1 


P 


28.7 ±0.9 


21.6 ±0.6 


13.2 ±0.4 


5.0 ±0.2 


0.73 ± 0.06 


P 


20.1 ± 1.0 


13.8 ±0.6 


9.2 ±0.4 


3.6 ±0.1 


0.47 ±0.05 



Pions dominate the charged particle multiplicity, but a large number of kaons and 
(anti)protons are produced. The inclusive yield of antiprotons is nearly comparable to 
that of protons. In the most central Au±Au collisions, the particle density at midrapidity 
(dN/dy) is ~ 20 for antiprotons and 28 for protons, not corrected for feed-down from strange 
baryons. 

The average transverse momenta increase with particle mass and with decreasing impact 
parameter. The mean transverse momentum increases with the number of participant nucle- 
ons by 20±5% for pions and protons, as shown in Figure El The (pt) of particles produced 
in p + p and pp collisions, extrapolated to RHIC energies, are consistent with the most pe- 
ripheral pion and kaon data; however, the {pt) of protons produced in Au±Au collisions is 
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TABLE III: The {px) in MeV/c for hadrons produced at midrapidity in each centrahty class. The 

errors are statistical only. The systematic uncertainties are 7%, 10%, and 8% for pions, kaons, and 
(anti) protons, respectively. 

0-5% 5-15% 15-30% 30-60% 60-92% 

7r+ 390 ± 10 380 ± 10 380 ± 20 360 ± 10 310 ± 30 

TT- 380 ±20 390 ± 10 380 ± 10 370 ± 20 320 ± 20 

K+ 560 ±40 580 ±40 570 ± 40 550 ± 40 470 ± 90 

K- 570 ± 50 590 ±40 610 ± 40 550 ± 50 460 ± 90 

p 880 ±40 870 ±30 850 ± 30 800 ± 30 710 ± 80 

p 900 ± 50 890 ±40 840 ± 40 820 ± 40 800 ± 100 



significantly higher. This dependence on the number of participant nucleons may be due to 
radial expansion. 
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FIG. 13: The integrated mean px for pions, kaons, and (anti)protons produced in the five different 

I ' 

classes of event centrahty 21]. The error bars are statistical only. The systematic uncertainties are 
7%, 10%, and 8% for pions, kaons, and (anti)protons, respectively. The open points are equivalent 
average transverse momenta from pp and pp data, interpolated to y/s=130 GeV. 
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C. TRANSVERSE MASS DISTRIBUTIONS 



Production of hadrons from a thermal source would make transverse mass the natural 
variable for analysis. Therefore we extract inverse slopes from the transverse mass distribu- 
tions by separately fitting a thermal distribution to each particle species. The Boltzmann 
distribution is given in equation 12. 



AmTe-"'^/^^ff, (12) 



2nmTd'mTdy 

We use a simple exponential, however, with no powers of in the prefactor, as shown in 
equation 10. This simplification is acceptable as the difference in the inverse slope is found 
to be less than 2%. The simple rriT exponential was also used in an equivalent analysis in 
Reference S^]. The inverse slope, Teff, can be compared to other experiments, provided the 
same momentum range of the spectrum is used for fitting. 

If the system develops collective motion, particles experience a velocity boost from this 
motion, resulting in an additional transverse kinetic energy component. This motivates use 
of the transverse kinetic energy, i.e. transverse mass minus the particle rest mass, for plotting 
data. Figure IT^ shows the transverse kinetic energy distributions (i.e. transverse mass minus 
the particle rest mass) for all positive particles (left) and negative particles (right). Pions 
are in the top panel, kaons in the middle panel, and (anti)protons in the bottom panel, with 
different symbols indicating different centrality bins. The solid lines are ttit exponential 
fits in the range (mj- — mg) < 1 GeV for all particle species while the dashed lines are the 
extrapolated fits. The pion spectra follow an exponential for 0.38 < (my — mo) < 1.0 GeV 
while the kaons and protons appear exponential over the entire measured m-p range. The 
same is true for the negative particles in the right panel; however, the antiprotons have more 
curvature for {rriT — mo) < 0.5 GeV. We extract T^ff by fitting exponentials of the form 
Equation El to the transverse mass spectra in the range (m^ — mo) <1 GeV. 

This range is chosen common for all particle species and minimizes contributions from 
hard processes. Caution must be taken when comparing T^ff values as the local slope of 
the transverse mass spectra varies somewhat over rriT for pions and antiprotons even within 
this fit range. The resulting values of T^ff for all particles and centralities are tabulated in 
Table US in units of MeV. The inverse slopes increase and then saturate for more central 
collisions for all particles except antiprotons. The fact that the inverse slope is different for 
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FIG. 14: Transverse mass distributions of pions (top), kaons (middle), and (anti)protons (bottom) 
in events with different centralities. Positive particle distributions are on the left and negative 
particle distributions are on the right. The solid lines are m-r exponential fits in the range (m-p — 
mo) < 1 GeV for all particle species while the dashed lines are the extrapolated fits. 

mesons and baryons and for central and peripheral events is consistent with the mean pt 
trends discussed above. 

We compare to published inverse slopes of transverse mass distributions at midrapidity 
from rriT exponential fits in the regiori (rnr — itlq) < 1^ GeV, listed in Table The 



comparison includes NA44 



27 



28 



29 



|30| and WA97 31, Q 
GeV; and, at y/s^N = 23 GeV at the ISR, Alper et al. [33[ and Guettler et al. |34i]. These 
data are chosen as they match the {rriT — rno) range used in fitting our data. For pions, 
the low-pT region of (m^ — mo) < 0.3 GeV, populated by decay of baryonic resonances, is 
systematically excluded from the fits. The effective temperatures are given in Table IVl with 
the references noted accordingly. 
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TABLE IV: The resulting inverse slopes in MeV after fitting an my exponential to the spectra in 
the range tjit — ttiq <1 GeV in each event centrality classes. The pion resonance region is excluded 
in the fits. The equivalent pT fit range for each particle is shown accordingly. The errors are 
statistical only. 

0-5% 5-15% 15-30% 30-60% 60-92% 

7r+ in 0.5 < pr < 1-05 GeV/c 216.8 ± 5.7 214.3 ± 4.6 217.4 ± 4.7 214.4 ± 5.2 176.9 ± 9.5 
TT- in 0.5 <pt < 1.05 GeV/c 215.8 ± 6.5 221.2 ± 5.6 225.3 ± 5.8 212.8 ± 5.7 215.8 ± 16.8 
K+ in 0.45 < pr < 1-35 GeV/c 233.2 ± 10.8 243.6 ± 9.8 242.4 ± 9.2 228.7 ± 10.2 182.3 ± 19.0 
in 0.45 <PT < 1.35 GeV/c 241.1 ± 15.8 244.5 ± 10.2 250.0 ± 12.3 224.2 ± 11.1 196.4 ± 22.3 
p in 0.55 <PT < 1.85 GeV/c 310.8 ± 14.8 311.0 ± 12.3 293.8 ± 11.4 265.3 ± 10.9 200.9 ± 14.8 
p in 0.55 <PT < 1.85 GeV/c 344.2 ± 25.3 344.0 ± 20.9 307.6 ± 17.1 275.1 ± 14.0 217.0 ± 28.3 



Radial flow imparts a radial velocity boost on top of the thermal distribution. Heavy 
particles are boosted to higher px, depleting the cross section at lower pt and yielding a 
higher inverse slope. Therefore, the observed inverse slope dependence on both centrality 
and particle mass implies more radial expansion in more central collisions. At CERN SPS, 
the Tgff depends on both mass and system size (the number of participating nucleons in the 
collision), indicating collective expansion. The Te// values at RHIC shown in Table IV are 
somewhat larger. 

In p-p collisions at similar y/s at the ISR, hadron spectra were analyzed in transverse 
mass, rriT, rather than transverse kinetic energy — mo To facilitate a direct 

comparison, figure El shows the PHENIX hadron spectra, including 7r° from the 10% most 
central Au+Au collisions. The spectra approach one another, but do not fall upon a universal 
curve, and thereby fail the usual definition of scaling. 

It has been suggested that at transverse mass significantly larger than the rest mass of the 
particle, thermal emission and radial flow may not be the only physics affecting the particle 
spectra. If heavy ion collisions can be described as collisions of two sheets of colored glass 
in which the gluon occupation number is sufficiently large to saturate, scaling of different 



hadron spectra with transverse mass is also predicted |23]- For Au+Au collisions at different 
impact parameters, the saturation scale differs, and some differences in the spectra may be 
expected. Nevertheless, the authors observe that the level of m-r scaling in our data is in 
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TABLE V: Inverse slope parameters (in MeV) of hadrons for p+p, p+nucleus, and central S+S, 
S+Pb, and Pb+Pb colliding systems at CERN energies. Pb+Pb is at ^^snn = 17 GeV, and the 
other systems at 23 GeV. The errors are statistical and systematic, respectively. 



Hadron 


Pb+Pb 


S+Pb S+S 


P+Pb 


p+S 


p+Be 


P+P 




156±6±23'' 


165+9+10'' 148+4+22'^ 


MS+S+IO** 


139+3+10'' 


148+3+10'' 


139+13+21*= 


K+ 


234±6±12'* 


181+8+10^ 180+8+9" 


172+9+10'' 


163+14+10'' 


154+8+10'' 


139+15+7" 


P 


289±7±14'^ 


256+4+10" 208+8+10" 


203+6+10'^ 


175+30+10" 


156+4+10" 


148+20+7" 


A 


289±8±29-^ 




203+9+205 








A 


287±13±29^ 




180+15+185 











"Reference 
^Reference 
"Reference 
'^Reference 
"Reference 
^Reference 
^Reference 



22 



(NA44 Collaboration). 
(T1A44 Collaboration). 
J (ISR). 

(NA44 Collaboration). 
(NA44 Collaboration). 
(WA97 Collaboration). 
(WA97 Collaboration). 



qualitative agreement with expectations from gluon saturation j37|]. Single particle spectra 
alone, however, are not sufficient to disentangle saturation from flow effects. 

It is often stated that rriT scaling holds in pp collisions at similar ^/s to RHIC [see data, for 
example, in references ^ and Scaling in uit, i.e. spectra following a universal curve 

in rriT, might be expected if the hadrons are emitted from a source in thermal equilibrium. 
It is instructive to note that reference 3^ states "Although the curves for different particles 
do come together, there is no real evidence for any universal behavior in this variable." Thus, 
scaling at the ISR was never claimed by the original authors. In central Au+Au collisions, 
the slopes and yields of vr, K and p approach each other as well, but figures IT^ and IT^ also 
do not support a truly universal behavior in m^. Therefore the apparent puzzle of how the 
data could exhibit both rriT scaling and the mass-dependent pt boost characteristic of radial 
fiow is no puzzle at all, as any "m^^ scaling" is only very approximate. 
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FIG. 15: The transverse mass distributions of inclusive identified hadrons produced in 10% central 
events, including the vr*^ as measured in the Electromagnetic Calorimeter in PHENIX and published 
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D. SUMMED CHARGED PARTICLE MULTIPLICITY 



As a consistency check we compare the measured rapidity densities as given in Section 
IIVBI to previously published pseudorapidity densities of charged particles. The measured 
dN/dy for each hadron species is converted to dN/di], and the total dN/dr/ is calculated 
by summation. Figure ^1 shows dN/d?7 per participant nucleon pair, compared to the 



measurement made by PHENIX usin g th e pad chambers alone [2^ as well as to PHOBOS 
and STAR yields in central collisions ■ We note that the lines correspond to the fit 

of a linear parameterization of Npart and Ncou to the PHENIX measurement (open circles) 
with a = 0.88 ± 0.28 and b = 0.34 ± 0.12 as described in 1411]. For the 5% central collisions, 
we measure 598 ±30, and is comparable to the STAR result of 567 ± 38 |40|, the PHOBOS 



result of 555 ±37 j42], and the PHENIX pad chamber result of 622 ±41 j26|. The agreement 
is excellent, allowing the results of this analysis to be used to decompose the particle type 
dependence of the charge particle multiplicity increase with centrality. 
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FIG. 16: The transverse mass distributions of the inclusive identified hadrons produced in the 5% 
central events. 

V. COMPARISON WITH MODELS 

A. HYDRODYNAMIC-INSPIRED FIT 



The charged particle pseudorapdity distributions are incompatible with a static thermal 



source, but the flat distribution observed in 



42| | reflects the strong longitudinal motion in the 



initial state. Consequently, the longitudinal momentum distribution is not an unambiguous 
sign of collective motion. Transverse momentum is, however, generated in the collision, so 
collective expansion may be more easily inferred from transverse momentum distributions. 

Following the arguments of the previous section, we analyze the particle rriT spectra. A 
parameterization of the rriT distribution of particles emitted from a hydrodynamic expanding 
hadron source is used. In order to determine the freeze-out temperature and collective flow 
without confusion from hard scattering processes, a limited px range is used in the fits. We 
include only particles with (m^ — mo) < 1 GeV in the fit. Pions with (m^ — rrio) < 0.38 
GeV are excluded to avoid resonance decays. All particles are assumed to decouple from 
the expanding hadron source j3| at the same freeze-out temperature, Tfo- This procedure 
allows us to extract Tjo and the magnitude of the collective boost in the transverse direction. 

The inverse slope includes the local temperature of a section of the hadronic matter along 
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FIG. 17: The transverse mass distributions of the inclusive identified hadrons produced in events 
with different centrahties. 

with its collective velocity. The simple exponential fit of Equation treats each particle 
spectrum as a static thermal source, and a collective expansion velocity cannot be extracted 
reliably from a single particle spectrum. However, the relative sensitivity to the temperature 
and collective radial flow velocity differs for different particles. By using the information 
from all the particles, the expansion velocity can be inferred. We fit all particle species 
simultaneously with a functional form for a boosted thermal source based on relativistic 
hydrodynamics j3] . 

Use of this form assumes that 

• all particles decouple kinematically on a freeze-out hypersurface at the same freeze-out 
temperature Tfo, 

• the particles collectively expand with a velocity profile increasing linearly with the 
radial position in the source (i.e., Hubble expansion where fluid elements do not pass 
through one another), and 

• the particle density distribution is independent of the radial position. 
Longitudinally boost invariant expansion of the particle source is also assumed. 
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FIG. 18: Both the total charged multiphcity (open) in References |39|, |40|, |4l| and the total iden- 
tified charged multiplicity (closed) scaled by the number of participant pairs are plotted together 
as a function of the number of participants. The lines correspond to the fit of a linear parameter- 
ization of Npart and Ncoii to the PHENIX charged multiplicity measurement (open circles) with 
a = 0.88 ± 0.28 and b = 0.34 ± 0.12 as described m^\. 

The transverse velocity profile is parameterized as: 

MO = /^r^r, (13) 

where ^ = and R is the maximum radius of the expanding source at freeze-out (0 < ,^ < 1) 
[4^ . The maximum surface velocity is given by and for a linear velocity profile, n = 

1. The average of the transverse velocity is equal to: 

/e^e "^^^ • ^ 

Each fluid element is locally thermalized and receives a transverse boost p that depends on 
the radial position as: 

p = tanh-i(/?T(0)- (15) 

The rriT dependence of the invariant yield ^^^^ is determined by integrating over the 
rapidity, azimuthal angle, and radial distribution of fluid elements in the source. This 
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procedure, discussed in Appendix O yields 

The parameters determined by fitting Equation 16 to the data are the freeze-out temperature 
Tfo, the normalization A, and the maximum surface velocity using a flat particle density 
distribution (i.e., /(^) = 1). 

To study the parameter correlations, we make a grid of combinations of temperature and 
velocity, and perform a chi-squared minimization to extract the normalization. A, for each 
particle type. The fit is done simultaneously for all particles in the range {mj- — uiq) < 1 
GeV. In addition to this upper limit in the fit, the pion fit range includes a lower limit of 
(mT — mo) > 0.38 GeV to avoid the resonance contribution to the low pt region (see Section 

The radial flow velocity and freeze-out temperature for all centralities are determined in 
the same way. The results are plotted together with the spectra in Figure El The hydro- 
dynamic form clearly describes the spectra better than the simple exponential in Figure IT^ 
The values for Tfo and (3^^°'^ are tabulated in Table IVIl 

TABLE VI: The minimum and the parameters Tfo and for each of the five centrality 

selections. The best fit parameters are determined by averaging all parameter pairs within the la 
contom'. The errors correspond to the standard deviation of the parameter pairs within the la \^ 
contom'. It is important to note that the fit range in Figure [T^ is the same as was used to fit 
exponentials to the spectra in Figure ITU 



Centrahty (%) 


xVdof 


Tfo (MeV) 


IJmax 


< /9t > 


0-5 


34.0/40 


121 ±4 


0.70 ±0.01 


0.47±0.01 


5-15 


34.7/40 


125 ±2 


0.69 ±0.01 


0.46±0.01 


15-30 


36.2/40 


134 ±2 


0.65 ±0.01 


0.43±0.01 


30-60 


68.9/40 


140 ±4 


0.58 ±0.01 


0.39±0.01 


60-92 


36.3/40 


161±11 


o.24±!]:i6 


0.16±[]i^ 



Figure 1201 shows contours for the temperature and velocity parameters for the 5% most 
central collisions. The n-sigma contours are labeled up to Sa. The contours indicate 
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FIG. 19: The parameterization and the px hadron spectra for all five centrality selections. 

strong anti-correlation of the two parameters. If the freeze-out temperature decreases, the 
flow velocity increases. The minimum is 34 and the total number of degrees of freedom 
(dof) is 40. The parameters that correspond to this minimum are Tfo = 121 ± 4MeV and 
pmax _ _|_ Q Qi 'jj^g quoted errors are the la contour widths of Ap^"-^ and ATfo. 
Within 3cT, the Tfo range is 106 - 141 MeV and the /^^'^^ range is 0.75 - 0.64. 

As a linear velocity profile (n = 1 in Equation IT^ is assumed, the mean flow velocity 
in the transverse plane is {/St) = 2/5™"^/3. If a different particle density distribution (for 
instance, a Gaussian function for f{^)) were used, then the average should be determined 
after weighting accordingly [44]. 

A similar analysis for Pb+Pb collisions at 158 A GeV, was reported by the NA49 Collab- 
oration in [4^ . Using the same hydrodynamic parameterization, simultaneous flts of several 
hadron species for the highest energy results in Tfo = 127 ± 1 MeV and Z?™""^ = 0.48 ± 0.01 
with x^lNDF = 120/43 for positive particles and Tfo = 114±2 MeV and p^""" = 0.50±0.01 
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with /NDF = 91/41 for negative particles (statistical errors only). Pions and deuterons 
are excluded from the fits to avoid dealing with resonance contributions to the pion yield 
and formation of deuterons by coalescence. The meson is included in the fit together with 
the negative particles. Previously, NA49 used a different parameterization to fit the charged 
hadron and deuteron spectra, as well as the rriT dependence of measured HBT source radii, 
resulting in overlapping contours with Tfo = 120 ± 12 MeV and = 0.55 ± 0.12 



1. VELOCITY AND PARTICLE DENSITY PROFILE 



In order to use and Tfo from the fits described above, one needs to know their 

sensitivity to the assumed velocity and particle density profiles in the emitting source. The 



choice of a linear velocity profi 
full hydrodynamic calculation 



e within the source is motivated by the profile observed in a 
1 

47| . which shows a nearly perfect linear increase of /5(r) with 
r. Nevertheless, we also used a parabolic profile to check the sensitivity of the results to 
details of the velocity profile. For a parabolic velocity profile (ra = 2 in Equation IT!^ . /?^"^ 
increases by ^ 13% and Tjo increases by ~ 5%. 

A Gaussian density profile used with a linear velocity profile increases jS^"-^ by ~ 2%, 
with a neglible difference in the temperature Tfo- As a test of the assumption that all the 
particles freeze out at a common temperature, the simultaneous fits were repeated without 
the kaons. The difference in Tfo is within the measured uncertainties. 



2. INFL UENCE OF RESONANCE PROD UCTION 



The functional forms given by Equations ^1 and El do not include particles arising from 
resonance or weak decays. As resonance decays are known to result in pions at low transverse 



momenta 



we place a pt threshold of 500 MeV/c on pions included in the 



hydrodynamic fit. A similar approach was followed by NA44, E814, and other experiments 
at lower energies, which performed in-depth studies of resonance decays feeding hadron 
spectra. However, these were for systems with higher baryon density, so we performed a 
cross check on possible systematic uncertainties arising from the pion threshold used in the 



fits. To estimate the effect of resonance decays were t 
calculate resonance contributions following Wiedemann 



ley not excluded from the fit, we 
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In order to reproduce the relative yields of different particle types, a chemical freeze-out 
temperature - different from the kinetic freeze-out temperature - and a baryonic chemical 
potential are introduced. Direct production and resonance contribution are calculated for 
pions and (anti)protons assuming a kinetic freeze-out temperature of 123 MeV, a transverse 
flow velocity of 0.612 (equivalent to {(3t) = 0.44), a baryon chemical potential of 37 MeV, 
and a chemical freeze-out temperature (when particle production stops) of 172 MeV. These 
parameters are chosen as they provide a reasonable description of the (anti)proton and pion 
spectra and yields (10% most central) and are in good agreement with chemical freeze- 
out analyses 0]. Most spectra from resonance decays show a steeper fall-off than the 
direct production, which should lead to a smaller apparent inverse slope, depending on what 
fraction of the low pt part of the spectrum is included in the fits. 

To measure the effect of resonance production on the spectral shape, the local slope is 
determined. For a given mrp bin number i, the local slope is defined as 

^ /-N rriTii + I) - rriTii - I) , . 

which is identical to the inverse slope independent of ttlt for an exponential. 
The difference in the local slope, 

\rp rpdirect rpincl 

local - J- local ' -I- locah [^°) 

is determined for direct and inclusive pions and (anti)protons. The differences are plotted 
as a function of rriT — rriQ in Figure |^ The difference in local slope for protons is below 13 
MeV for the full transverse mass range; the non-monotonic behavior for protons is caused 
by the relatively strong transverse flow. For pions, ATiocai decreases monotonically with rriT 
and is below 10 MeV above rriT = IGeV/c. A fit of an exponential to the pion spectra for 
{rriT — rriQ) > 0.38 GeV (which corresponds to pr > 0.5 GeV/c) yields a difference in inverse 
slope of 16 MeV with and without resonances. 



B. COMPARISON WITH HYDRODYNAMIC MODELS 

Hydrodynamic parameterizations as used in the previous Section rely upon many simpli- 
fying assumptions. Another approach to the study of collective flow is to compare the data 
to hydrodynamic models. Such models assume rapid equilibration in the collision and de- 
scribe the subsequent motion of the matter using the laws of hydrodynamics. Large pressure 
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FIG. 21: The difference in the local slope for direct and inclusive pions (solid) and (anti)protons 
(dashed). 

buildup is found, and we investigate this ansatz by checking the consistency of the data with 
calculations using a reasonable set of initial conditions. We compare to two separate models, 
the hydrodynamics model of Kolb and Heinz 5^ and the "Hydro to Hadrons" (H2H) 

model of Teaney and Shuryak [sgI . IstI i . The H2H model consists of a hydrodynamics calcu- 
lation, followed by a hadronic cascade after chemical freeze-out. The cascade step utilizes 
the Relativistic Quantum Molecular Dynamics (RQMD) model, developed for lower energy 
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heavy ion collisions 

In both models, initial conditions are tuned to reproduce the shape of the transverse 
momentum spectra measured in the most central collisions, along with the charged particle 
yield. Each model also includes the formation and decay of resonances. 

In the Kolb and Heinz model, the initial parameters are the entropy density, baryon num- 
ber density, the equilibrium time, and the freeze-out temperature which controls the duration 
of the expansion. The chemical freeze-out temperature is the temperature at which particle 
production ceases. The initial entropy or energy density and maximum temperature are 
fixed to match the measured multiplicity for the most central collisions using a parameteri- 
zation that is tuned to produce the measured dNch/(^'n dependence on both N^art and Ncou- 
A kinetic freeze-out temperature of Tjo = 128 MeV is used. Spectra from the Kolb-Heinz 
hydrodynamic model are shown in Figure |221 for pions (upper) and for protons (lower) as 
dotted lines. The solid lines are the results from the fits described in the previous sec- 
tions. The figure thus allows two comparisons. The similarity of the dashed and solid lines 
shows that the hydro dynamic-inspired parameterization used to fit the data results in a pr 
distribution similar to this hydrodynamic calculation. Comparing the dashed lines to the 
data points shows that the hydrodynamic model agrees quite well for most of the centrality 
ranges. It is important to note that the model parameters are uncertain at the level of 
10%, and, more importantly, the application of hydrodynamics to peripheral collisions may 
be less reasonable than for central collisions, as hydrodynamic calculations assume strong 



rescattering and a sufficiently large system size (discussed in 55|). 

In reference ^ |^ , the PHENIX p spectrum shape is well described by the H2H model 
with the LH8 equation of state. The cascade step in the H2H model removes the requirement 
that all particles freeze out at a common temperature. Thus the freeze-out temperature and 
its profile are predicted, rather than input parameters. Furthermore, following the hadronic 
interactions explicitly with RQMD removes the need to rescale the particle ratios at the 
end of the calculation, as they are fixed by the hadronic cross sections rather than at some 
particular freeze-out temperature. The LH8 equation of state includes a phase transition 



with a latent heat of 0.8 GeV. In j56|, |59[, the VL and the (f) are shown to decouple from the 
expanding system at T = 160 MeV, and they receive a fiow velocity boost of 0.45c. Pions 
and kaons decouple at T = 135 MeV with fiow velocity = 0.55c, while protons have T = 120 
MeV and fiow velocity > 0.6. These temperatures and fiow velocities are consistent with 
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FIG. 22: The hydrodynamics calculation with initial parameters tuned to match the most central 
spectra in the pT range 0.3 — 2.0 GeV/c. 

the values extracted from the data for the most central events. However, the average initial 
energy density exceeds the experimental estimate using formation time tq =1 fm/c. 

In Figure 1221 radial flow from the flts of the previous section are shown as a function 
of the number of participants for Tfo (top) and {/St) (bottom). There is a slight decrease 
of Tfo, while (Pt) increases with Npart, saturating at 0.45. The value of (Pt) from Kolb 
and Heinz is also shown, and agrees with the data reasonably well. In the plot of (/3t), the 
dashed line indicates the results of fltting the parameterization to the data while keeping 
Tfo flxed at 128 MeV to agree with the value used by Kolb and Heinz. Radial flow values for 
central collisions remain unchanged, while those in peripheral collisions increase. Even with 
the extreme assumption that all collisions freeze out at the same temperature, regardless of 
centrality, the trend in centrality dependence of the radial flow does not change. 
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FIG. 23: The expansion parameters Tjg and pjp'^^ as a function of the number of participants 
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As a comparison, the results from a hydrodynamic model calculation are also shown (Ref. 
55|). The dashed line corresponds to a fixed temperature of 128 MeV for all centralities in the 
parameterized fit to the data. 

C. HYDRODYNAMIC CONTRIBUTIONS AT HIGHER pr 

We use the parameters extracted from the fit to the charged hadron spectra in the low 
region to extrapolate the effect of the soft physics to higher px- This yields a prediction for 
the spectra of hadrons should a collective expanding thermal source be the only mechanism 
for particle production in heavy ion collisions. Comparing this prediction to the measured 
spectrum of charged particles or neutral pions should indicate the pt range over which soft 
thermal processes dominate the cross section. Where the data deviate from the hydrody- 
namic extrapolation, other contributions, as e.g. from hard processes or non-equilibrium 
production become visible. The approach described here differs from hydrodynamic fits to 



48 



I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I 1= 

CM 2h ^ pr=o.7o ^ 

= \ T,<, = 121 MeV 

> 

O 10 

O 

> 1 
■a 

I- -sT ^ 

Q.10 r o \\ 
= radial flow fit ji*+K*+p* 

CM 

-4 radial flow fit ji*+jt" 

^ I radial flow fit K%K" \ 

- radial flow f it p+ p 

.|Q I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I 

0.5 1 1.5 2 2.5 3 3.5 4 4.5 5 

(GeV/c) 




FIG. 24: The high pT hadron spectra in Reference |38| | compared to the fit results assuming radial 
flow from the vrKp spectra in the 5% central events. 



the entire hadron spectrum, as we fix the parameters from the low pt region alone, where 
soft physics should be dominant. 

The hadron spectrum is calculated using the fit parameters from the low pt region fits 
shown in the preceding section, and extrapolated to higher px- Figure |^ shows the calcu- 
lated spectrum for each particle type, and the sum of the extrapolated spectra is compared 
to the measured charged hadrons {h~^ + h~) in the 5% most central collisions. As non- 
identified charged hadrons are measured in t] rather than in y, the extrapolated spectra are 
converted to units of r/. This conversion is most important in the low p^ region. No addi- 
tional scale factor is applied - the extrapolation and data are compared absolutely. Below 
~ 2.5 GeV/c Pt, the agreement is very good, while at higher pt the data begin to exceed 
the hydrodynamic extrapolation. 

Other hydrodynamic calculations have been successful in describing the distributions 
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over the full pt range [60| with different parameter values. There are clear indications that 
particle production from a hydrodynamic source, if invoked to explain the spectra at low px, 
will have a non- negligible influence even at relatively large pt- Furthermore, the range of pt 
populated by hydrodynamically boosted hadrons is species dependent. This is clearly visible 
in Figure 1^ which shows that the extrapolated proton spectra have a flatter p^ distribution 
than the extrapolated pions and kaons. The yield of the "soft" protons reaches, and even 
exceeds, that of the extrapolated "soft" pions at 2 GeV/c pr- Therefore the transition from 
soft to hard processes must also be species dependent, and the boost of the protons causes 
the region where hard processes dominate the inclusive charged particle spectrum to be 
at significantly higher transverse momenta in central Au+Au than in p+p collisions. Our 
analysis suggests this occurs not lower than pj- = 3 GeV/c. 

D. HADRON YIELDS AS A FUNCTION OF CENTRALITY 

The previous discussion focused on the hadron spectra; now we turn to the centrality 
dependence of the pion, kaon, proton, and antiproton yields, which can shed further light 
on the importance of different mechanisms in particle production. It is instructive to see 
whether yields of the different hadrons scale with the number of participant nucleons, N parti 
the number of binary nucleon-nucleon collisions, Ncoiu oi some combination of the two. 

The total yields of the hadrons may be expected to be dominated by soft processes, and 
the wounded nucleon model of soft interactions suggests that the yields should scale as the 
number of participants, Npart- If each participant loses a certain fraction of its incoming 
energy, like e.g. in string models, where each pair of participants (or wounded nucleons) 
contributes a color flux tube, the total energy of the fireball formed at central rapidity would 
be proportional to the number of participants Npart- If, furthermore, the fireball is locally 
thermalized and particle production is determined at a single temperature, the multiplicity 
would scale with Npart- On the other hand, at very hig h p^, particle production may be 
dominated by hard processes and scale with A^^^;/|6ll. l62j|. 

In order to investigate the existence of scaling, the multiplicities are parameterized as: 

dN 

^ = C • (ATpart)"-- (19) 
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TABLE VII: Fit parameters for each particle species using equations 1191 and I2UI 



particle apart oi^ou 

7r+ 1.06 ±0.01 0.79 ±0.01 

TT" 1.08 ±0.01 0.80 ±0.01 

K+ 1.18 ±0.02 0.88 ±0.02 

K- 1.20 ±0.03 0.89 ±0.02 

J5+ 1.16 ±0.02 0.86 ±0.02 

p- 1.14 ±0.03 0.84 ±0.02 



and 

^ = C"-(iV,„,r^°"- (20) 



Fit results for these parameterizations are shown in Table rVTTl As can be seen, the exponents 
apart are > 1 for all species, while acou is consistently < 1. The production of all particles 
increases more strongly than with Npart, but not as strongly as with Ncoii- Small differences 
between the different particle species are apparent: The (anti-) proton yield increases more 
strongly than the pion yield, and the kaon yield shows the strongest centrality dependence. 
Remarkably, the yield fraction scaling beyond linear with Npart is larger for kaons, protons, 
and antiprotons than for pions. Perhaps it is not surprising that the yields do not scale 
simply with Npart', the collective flow seen in the pt spectra already shows that the nucleon- 
nucleon collisions cannot be independent. 

We next check whether the simple model of hadron yields can be brought into agreement 
with the data by adding a component of the yields scaling as the number of binary collisions, 
Ncoii- Such an admixture inspires simple two-component models j^,|6^. The nonlinearity 
of dN/dy on the number of participants is illustrated by the ratio (dN/dj)/ Npart, shown in 
Figure EHl as a function of centrality. The yields are seen to depend linearly on Ncou/Npart- 
As seen already from the exponents in Table IVTH the increase with centrality is strongest for 
kaons, intermediate for (anti-)protons, and weakest for pions. This indicates that protons 
and antiprotons have a larger component scaling with Ncoii than pions. 



We fit the yields per participant with Equation |^ As in 
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l62| we parameterize the 



multiphcity using two free parameters: Upp, the multiphcity in p±p collisions, and x, the 
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relative strength of the component scahng with N, 



coll • 



R = 



dN/dy 



L X } ' Tlpp — -\- X ' Tli^ 



part 



n, 



VP 



1 

2+^ 



coll 



N, 



part 



''pp 

2 



part 



(21) 



The results of the fit are shown as solid lines in Figure 031 The fit parameter values are given 
in Table IVlTTl All hadron species are well fit. The importance of the component scaling as 
Ncoii is largest for kaons and smallest for pions. 

TABLE VIII: Values of the parameters Upp and x from fitting Equation!^ to the observed dN/dy 

per Npart- 



n. 



pp 



IT ' 

K+ 
K- 

P 

P 



1.41 ±0.11 
1.10 ±0.11 
0.130 ±0.021 
0.089 ± 0.020 
0.089 ±0.013 
0.062 ±0.010 



0.028 ± 0.020 
0.085 ± 0.030 
0.232 ± 0.076 
0.326 ±0.132 
0.181 ±0.062 
0.172 ±0.068 



We check the consistency of the fits in Figure EHl with known hadron yields in p±p 
collisions by extrapolating the fits down to two participants (and one binary nucleon-nucleon 
collision). Isospin differences between p±p and Au±Au are ignored. The check is done by 
separately extrapolating the fitted fraction of yield which scales with Ncoii and the fraction 
scaling with Npart down to one nucleon-nucleon collision and two participant nucleons, and 
summing the result. One obtains particle ratios of K/tt = (8.7 ± 2.6)% and p/n = (4.9 ± 
0.8)%. These values fall between those measured at lower y/s at the ISR 0| and those at 
higher ^/s at the Tevatron [64], as expected since the RHIC energy lies in between. Thus 
the Au±Au data are shown to scale down to p±p reasonably. 

One may expect that the particle ratios at very high pt should be dominated by hard 
scattering, and therefore scale with the number of binary collisions. Consequently, we look 
at ratios of the Ncou scaling components alone, extrapolated down to one binary collision. 



The values are compared to measurements of hadron ratios at the ISR 



65| in Figures |2ni 
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FIG. 25: dN/ dy per participant of different particle species as a function of the number of collisions 
per participant. Kaon and (anti-)proton multiplicities are scaled by a factor of 20. 

and 123 The ratio of the extrapolated Au+Au yield fractions scaling as Ncou are shown as 
solid lines for > 2 GeV/c. The agreement with the p+p data at high pt is quite good. 

Finally, we directly compare p/vr and p/vr ratios in central Au+Au collisions with p+p, 
as a function of px- These ratios from the 10% most central data, using the charged particle 
measurement from this paper and neutral pions from 38|, are shown in Figure EHl The 
ratios show a steady increase up to 2.5 GeV/c in pt- Even though the simple extrapolation 
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FIG. 26: Kaon to pion ratio as a function of pT- The different points are measured in p+p collisions 



(data from ref. 



65|). The solid line is the asymptotic value for high pT in p+p derived from the 



hard scattering component of the fits using Equation to the measured centrality dependence 
of dN/dy in Au-Au collisions at -^/sjvtv = 130 GeV. The dashed lines indicate the corresponding 
uncertainty. 

of the Ncoii scaling yield fraction agreed with p+p, the ratios of the full yield significantly 
exceed those in the ISR measurements 6^. According to Gyulassy and collaborators, jGojl, 
this result may give insight into baryon number transport and the interplay between soft 
and hard processes. 

Of course, splitting the observed yields into portions which scale with Npart and Ncou 
is by no means a unique explanation of the data. The spectra and yields can also be 
well reproduced by thermal models, which break such simple scalings due to the multiple 
interactions suffered by the constituents. 

S.n,p,e .W.^ .ode. .hat .no. ._e and longitudinal flow H a. a.le to de- 
scribe the centrality dependence of the mid-rapidity vr^, K"^, p, and p yields by tuning the 
chemical freeze-out temperature Tch, the baryon chemical potential fiB and by introducing 
a strangeness saturation factor 7^. It was found that /i^ is independent of centrality, while 
both 7s and Tch increase from peripheral to central collisions. Within the same model, the 
centrality dependence of the particle yields at lower energy (a/s/vw" = 17 GeV 68, 0) are 
described by constant T^h and /i^. The strong centrality dependence in kaon production 
at both energies is accounted for by the increase in the strangeness saturation factor 7^. 
Although the integrated particle yields are very well described, such simple thermal models 
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FIG. 27: Antiproton to vr" ratio as a function of pT- The different points are measured va. p + p 
collisions (data from ref. 16^1 )• The solid line is the asymptotic value for high pT in p+p derived 
from the hard scattering component of the fits using Equation to the measured centrality 
dependence of dN/dy in Au-Au collisions at ^/snn = 130 GeV. The dashed lines indicate the 
corresponding uncertainty. 
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FIG. 28: a) Proton to pion ratio as a function oipT- b) Antiproton to pion ratio as a function 
of PT- The open circles represent measurements m. p + p colhsions (data from ref. y ). The filled 
circles show the 10 % most central Au+Au collisions. The neutral pion spectra are from the data 
pubhshed in (s^. 



do not attempt a comparison to the single particle spectra, which clearly indicate centrality 
dependent flow effects not included in the model. 

Thermal models which include hydrodynamical parameters on a freeze-out hypersurface 
to account for longitudinal and transverse flow can reproduce the absolutely normalized 
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particle spectra by introducing only two thermal parameters Tch and /is [70|, |7l|. In 
this approach, the thermal parameters are independent of centrality, while the geometric 
parameters are adjusted to reproduce the spectra. Good agreement with the data is obtained 
up to ~ 2 — 3 GeV/c, however an explicit comparison with the centrality dependence of 
the integrated mid-rapidity yields has not yet been made. 

This section shows that the yields of all hadrons increase more rapidly than linearly with 
the number of participants, but the increase is weaker than scaling with the number of binary 
collisions. The excess beyond linear scaling with Npart is strongest for kaons, intermediate 
for (anti-)protons, and weakest for pions. The centrality dependence of the total yields can 
be well fit with a sum of these two kinds of scaling. At high p^, the baryon and anti-baryon 
yields greatly exceed expectations from p+p collisions. Thermal models, which do not invoke 
strict scaling rules, can successfully reproduce the data as well, providing that they include 
the radial flow required by the px spectra. 



VI. SUMMARY AND CONCLUSION 



We have presented the spectra and yields of identified hadrons produced in ^/snn = 
130 GeV Au + Au collisions. The yields of pions increase approximately linearly with 
the number of participant nucleons, while the yield increase is faster than linear for kaons, 
protons, and antiprotons. 

Hydrodynamic analyses of the particle spectra are performed: the spectra are fit with a 
hydrodynamic-inspired parameterization to extract freeze-out temperature and radial fiow 
velocity of the particle source. The data are also compared to two full hydrodynamics 
calculations. The simultaneous fits of pion, kaon, proton, and antiproton spectra show that 
radial fiow in central collisions at RHIC exceeds that at lower energies and increases with 
centrality of the collision. The hydrodynamic models are consistent with the measured 
spectral shapes, extracted freeze-out temperature Tfo, and the fiow velocity Pt in central 
collisions. 

Extrapolating the fits to estimate thermal particle production at higher p^ allows us 
to study the soft-hard physics boundary by comparing to measured spectra at high pt- 
The yield of the "soft" protons reaches, and even exceeds, that of the extrapolated "soft" 
pions at 2 GeV/c pt- The sum of the extrapolated "soft" spectra agree with the measured 
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inclusive data to pt ~ 2.5 — 3 GeV/c. The transition from soft to hard processes must be 
species dependent, and the admixture of boosted nucleons imphes that hard processes do 
not dominate the inclusive charged particle spectra until approximately 3 GeV/c. 



APPENDIX A: DETERMINING Npart AND Ncoii 

As only the fraction of the total cross section is measured in both the ZDC and BBC 
detectors, a model-dependent calculation is used to map collision centrality to the number of 
participant nucleons, Npart, and the number of nucleons undergoing binary coUisions, Ncoii- 
A discussion of this calculation at RHIC can be found elsewhere |61|. 

Using a Glauber model combined with a simulation of the BBC and ZDC responses, Npart 
and Ncoii are determined in each centrality. The model provides the thickness of nuclear 
matter in the direct path of each oncoming nucleon, and uses the inelastic nucleon-nucleon 
cross section cr^^+^v determine whether or not a nucleon-nucleon collision occurs. We 
assume the following: 

• The nucleons travel in straight-line paths, parallel to the velocity of its respective 
nucleus. 

• An inelastic collision occurs if the relative distance between two nucleons is less than 

• Fluctuations are introduced by using the simulated detector response for both the 
ZDC and BBC. 



In this calculation, the Woods-Saxon nuclear density distribution (p(r)) is used for each 



nuc 
fm 



eus with three parameters, the nuclear radius R = 6.38Iqj3 fm, diffusivity d = 0.53±0.01 



3 



and the inelastic nucleon-nucleon cross section cr^^'^ — 40 ± 3 mb, 

Pir) = (Al) 
1 + 



APPENDIX B: INVARIANT CROSS SECTIONS 

Tabulated here are the measured invariant yields of pions, kaons, and (anti)protons pro- 
duced in Au+Au collisions at 130 GeV. The first set of tables (Tables lIXt IXIIjl are the 
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TABLE IX: Invariant yields for vr^, K^, and (anti)p measured in minimum bias events at midra- 
pidity and normalized to one rapidity unit. The errors are statistical only. 



(GcV/c) -± iv± p{->) 



0.25 


112±2 

109±2 








0.35 


56 ±1 
49.9 ±0.9 








0.45 


28.0 ±0.5 

24.1 ±0.5 




6.1±0.4 
4.6±0.4 




0.55 


15.7 ±0.3 




4.0±0.3 


2.3± 0.1 




14.6 ±0.3 




3.2±0.2 


1.2±0.1 


0.65 


9.1 ±0.2 




2.8±0.2 


1.8± 0.1 




8.7 ±0.2 




2.1 ±0.2 


1.17± 0.09 


0.75 


5.8 ±0.1 




1.7±0.1 


1.38± 0.08 




5.6 ±0.2 




1.6±0.1 


0.98±0.07 


0.85 


3.8 ±0.1 




1.30±0.08 


1.18±0.07 




3.6±0.1 




1.17±0.09 


0.95±0.07 


0.95 


2.40 ±0.08 




0.87±0.06 


0.98± 0.06 




2.28 ±0.08 




0.69±0.06 


0.65± 0.05 


1.05 


1.61 ±0.06 




0.62±0.04 


0.70± 0.04 




1.61 ±0.06 




0.53±0.05 


0.50± 0.04 


1.15 


1.03 ±0.04 




0.43±0.03 


0.60± 0.04 




1.17 ±0.05 




0.38±0.04 


0.35± 0.03 


1.25 


0.71±0.03 




0.33±0.03 


0.41± 0.03 




0.76 ±0.04 




0.27±0.03 


0.34±0.03 


1.35 


0.46 ±0.02 




0.20±0.02 


0.32±0.02 




0.54±0.03 




0.16±0.02 


0.22±0.02 


1.45 


0.35±0.02 




0.17±0.02 


0.23± 0.02 




0.31 ±0.02 




0.13±0.02 


0.18± 0.02 


1.55 


0.24 ±0.02 




0.10±0.01 


0.17± 0.02 




0.22 ±0.02 




0.10±0.01 


0.15± 0.02 


1.65 


0.16 ±0.01 
0.15 ±0.01 


58 


0.08±0.01 
0.07±0.01 





invariant cross sections plotted in Figures IHl and ^1 The second set of tables (Tables IXim 
IXVI|) are the invariant cross sections in equal pr bins as used in Figure ^2 



APPENDIX C: FREEZE-OUT SURFACE ASSUMPTIONS 



The freeze-out surface is a (r, 0, 77), where the radius r is between zero and R, the radius 
at freeze-out, the azimuthal angle is between zero and 27r, and the longitudinal space-time 
rapidity variable rj varies between —rjmax and rimax- In the Bjorken scenario, the freeze-out 
surface in space-time is hyperbolic, with contours of constant proper time r = ^/t"^ — z^. 
Assuming instantaneous freeze-out in the radial direction and longitudinal boost-invar iance, 
the model-dependence factors out of Equation [T2| and is included in the normalization 
constant A. 

At 130 GeV, the PHOBOS experiment measures the total charged particle pseudorapidity 
distribution to be flat over 2 units of pseudorapidity 4^ . The measured rapidity in PHOBOS 
is taken to be the same as the rapidity of the fireball, defined here as z. The rapidity variables 
in the integrand vanish for \z\ > 2. Therefore, the integration over the fireball rapidity is 
generally taken to be from —00 to +00 using the modified Ki Bessel function 

POO 

Ki{mT/T) = / cosh{z)e-'^'^'"''^^'^/'^dz (CI) 
Jo 

where the variable z is the fireball rapidity variable. The Ki bessel function can also result 
by integration over the measured rapidity y with the assumption that the freeze-out is 
instantaneous in the radial direction. In this case, no assumption is made on the shape 
of the freeze-out hypersurface. This also assumes that the total rapidity distribution is 
measured in the detector. What results is the single differential l/rriT dN/drnT-"^ 
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TABLE X: Pion invariant yields in each event centrality and pf bin measured at midrapidity, 
normalized to one rapidity unit. For each measured px bin, the positive pion cross section is the 
top row and the negative pion cross section is the bottom row. Errors are statistical only. 



(GeV/c) 


0-5% 


5-15% 


15-30% 


30-60% 


60-92% 


0.25 


355 ± 9 


282 ± 6 


186 ± 4 


81 ± 2 


13.2 ± 0.5 




371 ± 10 


275 ± 7 


180 ± 4 


74 ± 2 


12.1 ± 0.5 


0.35 


188 ± 5 


146 ± 3 


93 ± 2 


36.6 ± 0.8 


5.3 ± 0.2 




169 ± 5 


128 ± 3 


82 ± 2 


34.3 ± 0.9 


5.0 ± 0.2 


0.45 


95 ± 3 


74 ± 2 


48 ± 1 


17.5 ± 0.5 


2.7 ± 0.1 




86 ± 3 


63 ± 2 


40 ± 1 


15.7 ± 0.5 


2.1 ± 0.1 


0.55 


56 ± 2 


41 ± 1 


26.0 ± 0.7 


10.1 ± 0.3 


1.32 ± 0.09 




51 ± 2 


38 ± 1 


24.5 ± 0.8 


9.6 ± 0.3 


1.18 ± 0.09 


0.65 


32± 1 


25.6 ± 0.8 


15.0 ± 0.5 


5.3 ± 0.2 






30 ± 1 


22.6 ± 0.8 


14.5 ± 0.5 


5.7 ± 0.2 




0.70 










0.62± 0.04 
0.57± 0.04 


0.75 


21.1 ± 0.9 


15.4 ± 0.6 


9.9 ± 0.4 


3.6 ± 0.1 






20 ± 1 


15.3 ± 0.6 


9.5 ± 0.4 


3.5 ± 0.2 




0.85 


14.0 ± 0.7 


10.3 ± 0.4 


6.4 ± 0.3 


2.3 ± 0.1 






12.8 ± 0.8 


9.6 ± 0.5 


6.1 ± 0.3 


2.1 ± 0.1 




0.90 










0.19 ± 0.02 
0.24 ± 0.02 


1.00 


7.1 ± 0.3 


5.3± 0.2 


3.4 ± 0.1 


1.25 ± 0.05 






6.5 ± 0.3 


5.0 ± 0.2 


3.4 ± 0.1 


1.25 ± 0.05 




1.20 


3.2 ± 0.2 


2.2 ± 0.1 


1.47 ± 0.07 


0.55 ± 0.03 


0.064 ± 0.006 




3.3 ± 0.2 


2.6 ± 0.1 


1.51 ± 0.08 


0.63 ± 0.04 


0.061 ± 0.007 


1.40 


1.3 ± 0.1 


1.01 ± 0.07 


0.72 ± 0.04 


0.27 ± 0.02 






1.3 ± 0.1 


1.25 ± 0.09 


0.72 ± 0.05 


0.26 ± 0.02 




1.60 


0.55 ± 0.07 


0.57 ± 0.05 


0.33 ± 0.03 


0.14 ± 0.01 


0.015 ± 0.003 




0.51 ± 0.07 


0.57 ± 0.05 


0.30 ± 0.03 


0.12 ± 0.01 


0.014 ± 0.003 


1.80 


0.35 ± 0.05 


0.25 ± 0.03 


0.15 ± 0.02 


0.060 ± 0.008 






0.40 ± 0.06 


0.27 ± 0.03 


600.14 ± 0.02 


0.075 ± 0.009 




2.00 


0.18 ± 0.03 


0.13 ± 0.02 


0.10 ± 0.01 


0.023 ± 0.004 


0.005 ± 0.001 



TABLE XI: Kaon invariant yields in each event centrality and bin, measured at midrapidity 
and normalized to one rapidity unit. The top row in each pT bin is K^, and the bottom row is 
K~. Errors are statistical only. 



PT (GeV/c) 


0-5% 


5-15% 


15-30% 


30-60% 


60-92% 


0.44 










0.5 zb 0.1 
0.4 zb 0.1 


0.45 


21± 3 


16 ± 2 


10 zb 1 


4.3 zb 0.5 






21± 3 


13 ± 2 


7 zb 1 


2.5 zb 0.4 




0.54 










0.20 ± 0.07 

0.3 zb 0.1 


0.55 


15± 2 


11 zb 1 


6.6 zb 0.7 


2.4 zb 0.3 






13± 2 


8 lb 1 


4.8 zb 0.6 


2.3 zb 0.3 




0.65 


9 ± 1 


7.5 zb 0.7 


4.7 zb 0.4 


1.9 zb 0.2 






8 ± 1 


7.0 zb 0.8 


3.1 zb 0.4 


1.1 zb 0.2 




0.69 










0.18 zb 0.03 
0.14 zb 0.03 


0.75 


5.3 ± 0.7 


4.6 zb 0.5 


3.1 zb 0.3 


0.9 zb 0.1 






5.1 ± 0.8 


5.0 zb 0.6 


2.5 zb 0.3 


0.9 zb 0.1 




0.85 


5.7 ± 0.7 


3.6 zb 0.4 


2.1 zb 0.2 


0.67 zb 0.08 






3.6 zb 0.6 


3.7 zb 0.4 


2.0 zb 0.2 


0.64 zb 0.09 




0.89 










0.07 zb 0.02 
0.09 zb 0.02 


0.95 


3.0 ± 0.4 


2.3 zb 0.2 


1.5 zb 0.2 


0.51 ± 0.07 






2.4 ± 0.4 


2.3 zb 0.3 


1.0 zb 0.2 


0.35 zb 0.06 




1.05 


2.3 ± 0.3 


1.5 zb 0.2 


1.2 zb 0.1 


0.37 zb 0.05 






1.8 ± 0.3 


1.8 zb 0.2 


0.8 zb 0.1 


0.27 zb 0.05 




1.15 


1.6 lb 0.3 


1.3 zb 0.2 


0.62 zb 0.09 


0.29 zb 0.04 






1.4 ± 0.3 


0.9 ± 0.2 


0.7 ± 0.1 


0.23 zb 0.04 




1.17 










0.012zb 0.004 
0.015zb 0.005 


1.25 


1.1 ± 0.2 


1.0 zb 0.1 


0.66 zb 0.09 


0.15 zb 0.03 






1.2 ± 0.2 


0.8 ± 0.1 61 


0.41 ± 0.08 


0.15 ± 0.03 




1.35 


0.6 ± 0.1 


0.6 ± 0.1 


0.35 ± 0.05 


0.12 zb 0.02 





TABLE XII: (Anti)proton invariant yields in each event centrality and pT bin, measured at 
midrapidity and normalized to one rapidity unit. The top row in each pT is the proton cross 
section, and the bottom row the antiproton. The errors are statistical only. 



(GeV/c) 


0-5% 


5-15% 


15-30% 


30-60% 


60-92% 


0.545 










0.26 ± 0.06 
0.12 ± 0.05 


0.55 


8 ± 1 


4.9 ± 0.5 


4.0 ± 0.4 


1.6 ± 0.2 






4.2 ± 0.8 


2.8 ± 0.5 


2.0 ± 0.3 


0.8 ± 0.1 




0.65 


6.3 ± 0.7 


4.7 ± 0.4 


2.7 ± 0.2 


1.2 ± 0.1 






4.3 ± 0.7 


2.6 ± 0.4 


1.9 ± 0.3 


0.9 ± 0.1 




0.695 










0.14 lb 0.02 
0.14 ± 0.03 


0.75 


4.4 ± 0.5 


4.0 ± 0.4 


2.2 ± 0.2 


0.90 ± 0.08 






3.6 ± 0.5 


2.5 ± 0.3 


1.7 ± 0.2 


0.60 ± 0.08 




0.85 


3.9 ± 0.4 


3.3 ± 0.3 


1.9 ± 0.2 


0.75 ± 0.07 






2.9 ± 0.5 


2.5 ± 0.3 


1.8 ± 0.2 


0.62± 0.08 




0.895 










O.lOi 0.02 
0.06± 0.01 


1.00 


3.1 ± 0.2 


2.4 ± 0.2 


1.37 ± 0.09 


0.46 ± 0.03 






2.1 ± 0.2 


1.5 ± 0.1 


0.91 ± 0.08 


0.41 ± 0.04 




1.18 










0.031± 0.005 
0.018± 0.004 


1.20 


2.0 ± 0.2 


1.4 ± 0.1 


0.82 ± 0.06 


0.28 ± 0.02 






1.4± 0.2 


1.0 ± 0.1 


0.54 ± 0.05 


0.19 ± 0.02 




1.40 


1.1 ± 0.1 


0.74 ± 0.07 


0.46 ± 0.04 


0.14 ± 0.01 






0.9 ± 0.1 


0.50 ± 0.06 


0.32 ± 0.04 


0.13 ± 0.02 




1.58 










0.005± 0.002 
0.007± 0.002 


1.60 


0.54 ± 0.07 


0.49 ± 0.05 


0.25 ± 0.03 


0.09 ± 0.01 






0.49 ± 0.08 


0.37 ± 0.05 


0.20 ± 0.03 


0.053 ± 0.009 




1.80 


0.34 ± 0.05 


0.25 ± 0.03 


0.16 ± 0.02 


0.047 ± 0.007 





0.27 ± 0.06 0.11 ± 0.02 60.10 ± 0.02 



0.021 ± 0.005 



1.98 



0.003 ± 0.001 



TABLE XIII: Minimum bias invariant yields for all particles in equal pt bins. For each pt, the 
first line are the positive particle yields, and the second are the negative particle yields. The units 
are c'^/GeV^. 

PT{GeV/c) 7r± (anti)p 

0.25 112±2 

109±2 



0.35 


56 ±1 








49.9 ±0.9 






0.45 


28.0 ±0.5 


6.1±0.4 






24.1 ±0.5 


4.6±0.4 




0.55 


15.7 ±0.3 


4.0±0.3 


2.3±0.1 




14.6 ±0.3 


3.2±0.2 


0.38±0.02 


0.70 


7.3 ±0.1 


2.18±0.09 


1.55±0.06 




7.0 ±0.1 


1.9±0.1 


1.07±0.06 


0.90 


3.06 ±0.06 


1.07±0.05 


1.08±0.04 




2.89 ±0.07 


0.91±0.05 


0.79±0.04 


1.20 


0.91 ±0.02 


0.38±0.02 


0.49±0.02 




0.98 ±0.02 


0.32±0.02 


0.35±0.01 


1.60 


0.208 ±0.007 


0.104±0.006 


0.157±0.007 




0.193 ±0.007 


0.093±0.006 


0.119±0.007 


2.00 


0.050 ±0.003 




0.051±0.003 




0.053 ±0.003 




0.031±0.003 


2.45 


0.0028 ±0.0005 




0.013±0.001 




0.0034 ±0.0006 




0.009±0.001 


2.95 






0.0036±0.0006 








0.0022±0.0005 


3.55 






0.0007±0.0002 








0.0006±0.0002 
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TABLE XIV: Pion invariant yields in each event centrality normalized to one rapidity unit at 
midrapidity. The first line corresponds to positive pions, and the second to negative pions. 



Pt (GeV/c) 

If ± \ 1 / 


0-5% 


5-15% 


15-30% 


30-60% 


60-92% 


0.25 


355±9 


282±6 


186±4 


81±2 


13.2±0.5 




371±10 


275±7 


180± 4 


74 ± 2 


12.1 ± 0.5 


0.35 


188±5 


146ib3 


93± 2 


36.6± 0.8 


5.3 


± 0.2 




169±5 


128±3 


82± 2 


34.3± 0.9 


5.0 


± 0.2 


0.45 


95±3 


74±2 


48± 1 


17.5± 0.5 


2.7 


± 0.1 




86±3 


63±2 


40± 1 


15.7± 0.5 


2.1 


± 0.1 


0.55 


56±2 


41±1 


26.0± 0.7 


10.1± 0.3 


1.32 


± 0.09 




51±2 


38±1 


24.5± 0.8 


9.6 


± 0.3 


1.18 


± 0.09 


0.70 


26.3±0.8 


20.2±0.5 


12.3± 0.3 


4.4 


± 0.1 


0.62 


± 0.04 




24.7±0.9 


18.6±0.6 


11.8± 0.4 


4.5 


± 0.1 


0.57 


± 0.04 


0.90 


11.0±0.4 


7.9 


±0.2 


5.1 


± 0.2 


1.79 


± 0.06 


0.19 


± 0.02 




10.0±0.5 


7.9 


±0.3 


5.0 


± 0.2 


1.76 


± 0.07 


0.24 


± 0.02 


1.20 


3.1±0.1 


2.37 


±0.09 


1.50 


± 0.05 


0.58 


± 0.02 


0.064 


± 0.006 




3.4±0.2 


2.7 


±0.1 


1.67 


± 0.07 


0.64 


± 0.03 


0.061 


± 0.007 


1.60 


0.62±0.05 


0.54 


±0.03 


0.34 


± 0.02 


0.142 


± 0.009 


0.015 


± 0.003 




0.63±0.06 


0.58 


±0.04 


0.32 


± 0.02 


0.129 


± 0.009 


0.014 


± 0.003 


2.00 


0.17±0.02 


0.14 


±0.01 


0.083 


± 0.009 


0.027 


± 0.003 


0.005 


± 0.001 




0.20±0.03 


0.14 


±0.02 


0.09 


± 0.01 


0.035 


± 0.004 


0.004 


± 0.001 


2.45 


0.005 ±0.003 


0.009 


±0.003 


0.006 


± 0.002 


0.0017 


± 0.0006 








0.011±0.005 


0.011 


±0.004 


0.005 


± 0.002 


0.0025 


± 0.0009 







College of Arts and Sciences, Vanderbilt University (U.S. A), Ministry of Education, Culture, 
Sports, Science, and Technology and the Japan Society for the Promotion of Science (Japan), 
Russian Academy of Science, Ministry of Atomic Energy of Russian Federation, Ministry 
of Industry, Science, and Technologies of Russian Federation (Russia), Bundesministerium 
fuer Bildung und Forschung, Deutscher Akademischer Auslandsdienst, and Alexander von 
Humboldt Stiftung (Germany), VR and the Wallenberg Foundation (Sweden), MIST and 
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TABLE XV: Kaon invariant yields in each event centrality normalized to one rapidity unit at 
midrapidity. The first line corresponds to positive kaons, and the second to negative kaons. 



Pt (GeV/c) 


0-5% 


5-15% 


15-30% 


30-60% 


60-92% 


0.45 


21±3 


16± 2 


lOib 1 


4.3 ± 0.5 


0.5 


± 0.1 




21±3 


13± 2 


7 ± 1 


2.5 ± 0.4 


0.4 


± 0.1 


0.55 


15±2 


llzb 1 


6.6 ± 0.7 


2.4 ± 0.3 


0.20 


± 0.07 




13ib2 


8 ± 1 


4.8 ± 0.6 


2.3 ± 0.3 


0.3 


± 0.1 


0.70 


8.0±0.7 


6.6 ± 0.5 


4.3 ± 0.3 


1.5 ± 0.1 


0.18 


lb 0.03 




7.0±0.8 


6.5 ± 0.6 


3.1 ± 0.3 


1.1 ± 0.1 


0.14 


± 0.03 


0.90 


4.5±0.4 


3.1 ± 0.2 


1.9 ± 0.1 


0.62 ± 0.06 


0.06 


± 0.02 




3.3±0.4 


3.3 ± 0.3 


1.6 ± 0.2 


0.5 ± 0.06 


0.09 


± 0.02 


1.20 


1.4±0.1 


1.10 ± 0.08 


0.68 ± 0.05 


0.22 ± 0.02 


0.012 


± 0.004 




1.3±0.1 


0.97 ± 0.08 


0.50 ± 0.05 


0.17 ± 0.02 


0.015 


± 0.005 


1.60 


0.36±0.05 


0.29 ± 0.03 


0.17 ± 0.02 


0.062 ± 0.007 


0.008 


± 0.002 




0.29±0.05 


0.27 ± 0.03 


0.20 ± 0.02 


0.058 ± 0.008 


0.004 


± 0.002 



the Natural Sciences and Engineering Research Council (Canada), Conselho Nacional de 
Desenvolvimento Cientifico e Tecnologico and Fundagao de Amparo a Pesquisa do Estado 
de Sao Paulo (Brazil), Natural Science Foundation of China (People's Repubhc of China), 
Centre National de la Recherche Scientifique, Commissariat a I'Energie Atomique, Insti- 
tut National de Physique Nuclaire et dc Physique des Particules, and Association pour la 
Recherche et le Developpemcnt dcs Methodes et Processus Industrials (France), Department 
of Atomic Energy and Department of Science and Technology (India) , Israel Science Foun- 
dation (Israel), Korea Research Foundation and Center for High Energy Physics (Korea), 
the U.S. Civihan Research and Development Foundation for the Independent States of the 
Former Soviet Union, and the US-Israel Binational Science Foundation. 
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TABLE XVI: (Anti)proton invariant yields in each event centrality normalized to one rapidity 
unit at midrapidity. The first line corresponds to protons, and the second to antiprotons. 



(Gc'V/c) 


0-5% 


5-lo 




15-31 


1% 


30-60% 


GO-92'X 


0.55 


8±1 


4.9 ± 0.5 


4.0 ± 


0.4 


1.6 ± 


0.2 


0.26 ± 0.06 




4.2±0.8 


2.8± 


0.5 


2.0± 


0.3 


0.8± 


0.1 


0.12± 0.05 


0.70 


5.4±0.4 


4.5± 


0.3 


2.5± 


0.2 


1.06± 


0.07 


0.14± 0.02 




4.4±0.5 


2.9± 


0.3 


2.0± 


0.2 


O.SOib 


0.08 


0.14± 0.03 


0.90 


3.9±0.3 


3.1± 


0.2 


1.9± 


0.1 


O.Tlib 


0.05 


O.lOi 0.02 




2.8±0.3 


2.1± 


0.2 


1.4± 


0.1 


0.58± 


0.05 


0.06± 0.01 


1.20 


l.QibO.l 


1.37± 


0.08 


0.78± 


0.04 


0.26± 


0.02 


0.031± 0.005 




1.3±0.1 


0.96± 


0.07 


0.56± 


0.04 


0.21± 


0.02 


0.018± 0.004 


1.60 


0.60±0.06 


0.44± 


0.03 


0.27± 


0.02 


0.087± 


0.008 


0.005± 0.002 




0.49±0.06 


0.34± 


0.03 


0.19± 


0.02 


0.062ib 


0.007 


0.007± 0.002 


2.00 


0.20±0.03 


0.15± 


0.02 


0.09± 


0.01 


0.025± 


0.003 


0.003± 0.001 




0.15±0.03 


0.07± 


0.01 


0.055± 


0.009 


0.019 ± 


0.003 


0.0005± 0.0005 


2.45 


0.06±0.01 


0.040± 


0.007 


0.020 ± 


0.004 


0.006± 


0.001 


O.OOlOib 0.0006 




0.04±0.01 


0.028± 


0.006 


O.Olli 


0.003 


0.005± 


0.001 


0.0020 ± 0.0008 


2.95 


0.015±0.005 


O.OOTi 


0.002 


0.005± 


0.002 


0.0023± 


0.0007 


0.0006± 0.0004 




0.013±0.005 


0.008± 


0.003 


0.002± 


0.001 


0.0003± 


0.0003 


0.0005± 0.0003 


3.55 


0.003 ±0.002 


0.0012 ± 0.0007 


0.0014± 


0.0006 


0.0005± 


0.0002 






0.002±0.001 


O.OOllzb 0.0007 


0.0013± 


0.0006 


0.0002± 


0.0002 
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